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Introduction 


Ball bearing fatigue data consist of a list of the number of hours endured by the individual 
bearings of a group of a certain size. We may be given, for example, the number of hours 
that each of ten identical bearings endured under a fixed load. From this sample of ten bear- 
ings we would like to determine the general distribution of the fatigue life of such bearings. 
In other words, we would like to determine the nature of the population, given the sample. By 
the population we mean the large general group of similar bearings of which the sample is 
only a small portion. Thus, the population may consist of all the bearings of the same type 
as those of the sample which were produced during a certain month, or in a certain plant, or 
by a special process. 


The usual method of estimating populations from samples is to put the data into a frequency 
table and then to construct a histogram, or frequency polygon, which is finally smoothed into 
a frequency curve. Such a procedure, however, is not reliable for small samples, since the 
results vary too much with a change in the size of the class interval, and generally, false 
peaks cannot be distinguished from true ones. With the small samples which arise in fatigue 
life studies it has been found to be much more fruitful to order the data according to life and 
to construct a cumulative plot, or ogive curve. In other words, the estimation of the integral 
of the frequency function is easier than the direct estimation of the frequency function itself. 


In order to make a cumulative plot, it becomes necessary to decide what rank is to be as- 
Signed to each particular failure in a group. The first failure in a group of ten tested will 
have a definite percentage of the population failing before it, if the entire population were 
tested. If we knew exactly the percentage of the population below the first failure in ten, then 
that percentage would be the true rank of the first failure in ten. However, since we do not 
know the true rank, the best that we can do is to estimate what it is. We use an estimate such 
that in the long run the positive and negative errors of the estimate cancel each other. That 
is, half of the time we would give the first failure in ten a rank that is too high, and the other 
half of the time we would give it a rank that is too low. A rank with this property is called a 
“median” rank. 


It is the purpose of the present mathematical discussion to determine the median ranks to be 


assigned respectively to the first, second, third, etc., up to the tenth failure in a sample of 
ten failures, and, likewise, for all the failures in a sample of any finite size. 


Deriving the Median Ranks 





Assume the following situation is given: 
A sample of n observations in numerical order: 0%» o%» o%s +e » oXn: 
Probability density function of population: f(x) (unknown) 


Cumulative distribution function of population: F(x) (unknown) 
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We define, 


True rank of ,X, 


, = F(,X,) (unknown) 


Z 
True rank of a oe 


2 = F(,X,) (unknown) 


True rank of ,X, = ,Z, = F(,X,) (unknown) 


In general, True rank of oX, = nb; = P(oX; ) (unknown) 


Since the true rank of an observation is unknown, the best that we can do is to estimate what 


that true rank is. Consider next the set of all possible samples of size n from the same 
population. The table below is a partial list of these samples: 


(1) » sae _ aa 
Memeges Bs ake pete polkas’ °° * * * *vem 
~ — = en y ‘?) 
Sample $: XI. eke seke 5° * °° * * *yeke 
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Sample r: ome - * alle oo & a: * sale 


The true ranks of the jth observations in these samples of size n may be listed as follows: 
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We shall show that these ranks are distributed according to the probability density function 


n Sai 


: 9 1 .2 
(j-1)! (n-j)t 79 (1-n 45 ) 


é(nZ;) = "$ 





This follows from the multinomial theorem, as demonstrated below: 


Divide the entire population into the following three regions: 


Regio 





The expressions F (Xj), {(ox;) dx, and 1 - F (Xj) represent the probabilities that a single ob- 


Servation would fall into regions 1, 2, and 3, respectively. 


Therefore, it follows that the probability of exactly (j-1) values falling into region 1, and 
exactly one value falling into region 2, and exactly (n-j) values falling into region 3 is given 
by the expression 


n! 


j-1 ‘ -j 
Gay ater (F(X, )]°~*+ £(gX, Jax * [1 - F(,X,)]” 


(This is an application of the multinomial theorem on joint probabilities.) 


Now put F(,X;) = nd; , Then, f(,X, )dx = qd, Z, 


Hence, the above probability becomes 








CUP Gacy M9 (Fyn 
Therefore, the probability density function of n2j is 
E602) = ap targy Ms = 02)" (A) 
By definition, the mean value of nZj; is 
1 1 
we °* arc C=) es Be 
i) ° 


This integral can be reduced to a Beta function, and its value turns out to be j/(nt+1). This is 
the mean rank of the jth observation in n. We could use this as an estimate of the true rank 
of the jth observation in n, but in certain applications the median value of nZj has been found 
to be more useful an an estimate. Let us now proceed to find the median value of nZj; i.e., 
the median rank of the jth value in n. 


By integration of equation (A), we find that the cumulative distribution function of nZj is 





(2) ( j-1) 
a/ se. fs a + py \ael n 2 gy \a~2 oF n j-1 n-j+i 
G(, 25) 1 (1 ny) My Z; (1-2, ) “= nZ,(1-n2,) oo "Oca 44 (1-,,2, ) (B) 
(Note: es = n(n-1) a’? = n(n-1)(n-2), etc.,.. ,?°") = n(n-1)(n-2) ********(n-j+2).) 
The median value of nZj is found by putting G(Z;) = in equation (B) and then solving for n2jp 
i.e., by solving the following equation for nZj: 
(2) (j-1) 

nei n_ ,t he elite CLE . ginmleg 7 yerstt 

1 = (1-n 25)" -tn25(1-n25) = nZy(Im-w4y) 0 =" = (j-1)) "3 (1-n25) 3 (C) 


Such an equation contains exactly one real root between 0 and 1, and Table 1 consists of such 
roots for n and j up to 20. The table is constructed by solving equation (C) exactly for j = 1 
and for j = n, for each value of n up to 20. The intermediate values in each column are then 
filled in by forming an arithmetic progression. When this is done, we find that in every case 
495 *G(,,2)) = .505, which is good enough for practical purposes. 








For n>20 the following convenient formula may be used: 


j - (1-In 2) = (21a 2 -1) (>) 
Aa (j) = , where ) ,,(j) denotes the median rank of the jth 





value inn. 


An Application 


Eighteen identical ball bearings were placed in test machines and subjected to a fixed radial 
load. The following table lists the number of hours the individual bearings endured at 2000 
rpm: 





TABLE A 
Failure No. Endurance in Hours 
1 183 
2 355 
3 538 
4 618 
5 697 
6 834 
7 862 
8 887 
i) 1056 
10 1147 
11 1351 
12 1506 
13 1578 
14 1607 
15 1683 
i6 1710 
17 2020 
18 2410 


As a result of testing thousands of bearings, it has been found that the fatigue life follows a 
Weibull distribution, whose probability density function is 


B-1 B :, 
f(x) = aBx exp(-ax ) (a and B are positive constants.) 
Therefore, the cumulative distribution function of the life is 
p 
F(X) = 1- exp(-ax ), where F(x) is the fraction failed within the first x hours. 


Let S(x) = the fraction surviving the first x hours 


Then, S(x) = exp(-ax 


= oxp(ax” 
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ln tne = Ina +Blnx 





Hence, if we plot Log at vs. x on logarithmic paper, we should get a straight line, if the 


distribution is a Weibull distribution. 
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Figure l, 


Figure 1 shows a plot of the data in Table A. It is plotted on some special probability paper 
on which a Weibull distribution becomes a straight line of slope B. The vertical scale repre- 


sents Log Log on = Log Log iE but is graduated in terms of values of F(x), the frac- 


tion failed within the first x hours, while the horizontal is an ordinary logarithmic scale, 
representing Log x. 


In the left half of the figure, each individual bearing life is given its mean rank, while in the 
right half of the figure, it is given its median rank. The same straight line is drawn in both 
cases, and it represents the maximum likelihood estimate of the Weibull population from 
which the sample is drawn. Determining the maximum likelihood estimate from a set of data 
as shown in Table A is an arduous task, and it would be desirable to have a quick graphical 
method of determining the best position of the straight line representing the Weibull population. 


To draw the line very near the lower extreme values when mean ranks are used would lead to 
a slope which is too small, because of the high probability of lower extreme values falling 
considerably to the left in such a case. On the other hand, it has been found that if median 
ranks are used in the plotting, then the danger of underestimating the slope is eliminated be- 
cause, in that case, a point is just as liable to fall to the right as to the left of the maximum 
likelihood line. In other words, we can very quickly arrive at an estimate of the population 
parameters by drawing a line by sight which takes the general direction of the array of points 
and which splits the array 50-50. Forthis reason median ranks are preferred to mean ranks. 


Additional Remarks 





Equation (B) can also be used to determine tolerance limits for the jth value inn. If we were 
seeking the rank below which we can expect 30% of all jth observations in n, we would solve 


the following equation for nZj : 


(2) (j-1) 
— n-1 n . n-2 j-1 
1 -(1-, 23) Dy Zs (1-nZ,) - rT n Zs (1-n2;) —***= Ti) nZ, (1-, 25) 10 











Similarly, if we were seeking the rank below which we can expect 70% of all jth observations 
in n, we would solve the following equation for ,Z;: 








ji’ 
4 ( 7 )* 7 (1 7 , 83 a” 2° Z a-2 eece ” pee z?~* (4. Z, yprern oa (E) 
. ~(1-n 4, Ry “, “ni “4 = 2! n y (lon 3? - - G-1)!" 3 nj 10 


Figure 2 illustrates the application of this process to the data in Figure 1. Solutions of equa- 
tions (D) and (E) for n = 18 are found in the columns headed by 18 in Tables 2 and 3 respec- 
tively. 
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Figure 2, 

Referring to Tables 2 and 3 for n = 18 and j = 2, we see that they tell us that between the 
ranks .0608 and .1302 we can expect the middle 40% of all second failures in eighteen. On the 
straight line in Figure 2, 280 hours corresponds to the rank .0608 and 440 hours corresponds 
to the rank .1302. Hence, the middle 40% of all second failures in eighteen will lie between 
280 hours and 440 hours, if such samples of eighteen come from the population represented 
by the straight line. The same type of limits can be determined for the first failure, the 
third failure, the fourth failure, etc., up to the eighteenth failure in the sample of eighteen. 
When this is done, we obtain a middle 40% band as shown in the right half of Figure 2. 


A complete list of solutions of equation (D) for n and j up to 20 is given in Table 2. For n>20, 
the following formula may be used to find the 30% rank of the jt value in n: 


i. 
" n+1 


125 M-M) nt] 2(1-2M), / nt2 _ 
ck oo _ A ava ” A 
pals) © B ~{.068 sts iE he Ss Se n+3 = 











The above formula applies only to the lower half of a sample. For the upper half of the 
sample, use the following formula: 





4 M(1-M) n+1 
=)") 4 c 


aR (ntl-j) = 1-M -(1.01- .49( = 


(M and a, are the same as above.) 
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Likewise, a complete list of solutions of equation (E) for n and j up to 20 is given in Table 3. 
For n>20, the following formulas may be used for the 70% rank of the jth value in n: 











For the lower half of a sample: Rn (J) 2 j- 3 Rn (nt1-)) = M +[1.01- 49 (= fe A a 
° e nt 
For the upper half of a sample: ,R,(nt1-j) = 1- 3Rn (3) = 1 -M +(.645 - 125) M(1-M) 
‘ 5s nt+2 
+7 
(In all these formulas, J < =, and M and a, are the same as previously defined.) 
TABLE 1 — MEDIAN RANKS 
Sample Size = n 

A 1 2 3 4 5 _ 6 7 8 9 10 
1 .5000 .2929 .2063 .1591 .1294 .1091 .0943 .0830 .0741 .0670 
2 .7071 .5000 .3864 3147 .2655 .2295 .2021 .1806 .1632 
3 .7937 .6136 .5000 4218 .3648 .3213 .2871 .2594 
4 .8409 .6853 .5782 .5000 .4404 .3935 .3557 
5 .8706 .7345 .6352 .5596 .5000 4519 

1 6 8909 .7705 .6787 .6065 .5481 
7 .9057 .1979 .7129 .6443 
8 .9170 .8194 .7406 
9 .9259 .8368 
10 .9330 

Sample Size = n 

4 11 12 13 14 _15 16 17 18 19 20 
1 .0611 .0561 .0519 .0483 .0452 .0424 .0400 .0378 .0358 .0341 
2 1489 —_.1368 .1266 .1178 1101 .1034 .0975 .0922 .0874 .0831 
3 .2366 .2175 .2013 .1873 1751 .1644 .1550 .1465 .1390 .1322 
4 .3244 .2982 .2760 .2568 .2401 .2254 .2125 .2009 .1905 .1812 
5 4122 .3789 .3506 .3263 .3051 .2865 .2700 .2553 .2421 .2302 
6 .5000 4596 4253 .3958 .3700 .3475 .3275 .3097 .2937 .2793 
7 .5878 .5404 .5000 .4653 .4350 4085 .3850 .3641 .3453 .3283 
8 .6756 .6211 5747 5347 .5000 4695 4425 .4184 .3968 .3774 


9 .7634 .7018 .6494 .6042 .5650 5305 .5000 4728 4484 4264 
. 10 8511 .7825 .7240 6737 .6300 5915 5575 5272 .5000 4755 
yu .9389 .8632 .1987 -7432 .6949 .6525 .6150 .5816 .5516 5245 


12 .9439 .8734 8127 .1599 .7135 .6725 .6359 .6032 .5736 
13 .9481 .8822 .8249 .1746 .7300 .6903 .6547 .6226 
14 9517 .8899 .8356 .1875 “1447 .7063 6717 
15 .9548 .8966 .8450 .7991 .1579 .1207 
16 .9576 .9025 .8535 .8095 .7698 
17 .9600 .9078 .8610 .8188 
18 .9622 .9126 .8678 
19 .9642 .9169 
20 .9659 
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TABLE 2 — 30% RANKS 


Sample Size 


=n 











1 2 3 os 5 6 7 8 9 10 
~ .3000 .1633 1121 .0853 .0689 .0577 .0497 .0436 .0389 .0350 
5477 .3633 .2724 .2180 .1818 .1559 .1365 .1214 .1093 
.6694 4916 .3898 .3233 .2763 .2413 .2142 .1926 
.7401 .5780 4761 .4053 .3533 .3128 .2802 
.7860 .6396 5414 A701 4157 .3726 
.8182 .6857 .5925 5225 4675 
.8400 .7214 .6335 .5655 
.8603 .7499 .6670 
.8748 -TT31 
.8866 

Sample Size =n 

11 12 13 14 15 16 17 18 19 20 
~ 0319 ~—«—«.0293 .0271 .0252 .0235 .0220 .0208 .0196 .0186 0177 
.0994 .0911 .0841 .0781 .0729 .0684 .0644 .0608 .0576 .0547 
.1749 .1600 .1480 .1370 .1280 .1200 .1130 .1066 .1010 .0959 
.2555 .2336 .2154 .1987 .1854 1737 .1635 .1543 .1462 .1389 
.3392 .3097 .2848 .2626 .2447 .2292 .2156 .2035 .1928 .1834 
4249 .3878 .3560 .3284 .3056 .2861 .2689 .2540 .2405 .2291 
5125 ABTT 4289 .3957 .3679 .3442 .3234 .3050 .2891 .2757 
.6021 .5494 .5035 4645 .4316 .4035 .3789 .3572 .3385 .3229 
.6949 .6329 .5798 5348 4967 .4640 .4354 .4100 .3886 .3706 
.7923 .7186 .6581 .6066 .5632 5257 .4929 4642 4394 .4188 
.8963 .8086 .7388 .6804 .6311 .5886 .5514 .5182 .4909 4675 
.9045 .8225 -7564 .7006 .6527 .6109 .5748 .5431 .5167 
9115 .8345 -TT17 .7183 .6716 .6307 .5960 .5664 
.9176 .8450 .7853 .7335 .6888 .6497 .6166 
.9229 .8543 .1973 -7473 .7043 .6674 
.9275 .8624 .8080 -7600 .7191 
.9316 .8698 8177 .7720 
9353 .8764 .8264 
.9386 .8824 
.9416 





TABLE 3 — 70% 


RANKS 


Sample Size = n 











1 2 3 4 5 6 1 8 g 10 
1 .7000 .4523 3306 .2599 .2140 .1818  .1600  .1397 .1252 .1134 
2 .8367 .6367 .5084 4220 .3604 .3143 .2786 2501 .2269 
3 8879  .7276 .6102 .5239 .4586  .4075 .3665  .3330 
4 9147 .7820 .6767 .5947 .5299 .4775  .4345 
5 9311 .8182 .7237 .6467 .5843 .5325 
6 9423 .8441 .7587 .6872 .6274 
1 9503 .8635  .7858 .7198 
8 9564 .8786 .8074 
9 9611 .8907 

10 .9650 
Sample Size = n 
11 12 13 14 15 16 17 18 19 20 

~ | (1037 .0955 .0885 .0824 .0771 .0725 .0684 .0647 .0614 # .0584 
2 .2077 .1914 .1775 .1655 .1550 1457 .1376 .1302 .1236 © .1176 
3 .3051 .2814 .2612 .2436 .2283 2147 .2027 4.1920 .1823 .1736 
4 .3979 .3671 .3419 3196 .2994 4.2817 .2665  .2527 .2400 .2280 
5 .4875 .4506 .4202 .3934 .3689 .3473 .3284 #£+4.3112 .2957 .2809 
6 .5751 .5323 .4965  .4652 .4368 4114 .3891 .3693 .3503 .3326 
7 (6608 .6122 .5711 .5355 .5033 4743 .4486 .4252 .4040 .3834 
6 .7445  .6903 .6440 .6043 .5684 5360 .5071 .4818  .4569 .4336 
9 .8251 .7664 .7152 .6716 .6321 5965 .5646 .5358  .5091 .4833 

10 .9006 .8400 .7846 .7374 4.6944 6558 4.6211 .5900 .5606 .5325 
11 .9681 .9089 .8520 .8013 .7553 7139 .6766 .6428 4.6114 &# .5812 
12 9707 .9159 .8630 .8146 .7708 .7311 .6950 .6615 .6294 
13 9729 .9219 .8720 .8263 .7844 .7460 .7109 .6771 
14 9748  .9271 .8800 .8365 .7965  .7595  .7243 
15 9765  .9316 .8870 .8457 .8072 .7709 
16 9780 .9356 .8934 .8538 .8166 
17 9792 .9392  .8990 .8611 
18 9804 .9424 .9041 
19 9814 .9453 
20 .9823 








APPLICATIONS OF ELEMENTARY MATHEMATICS 
TO PROBLEMS IN MAGNETISM 


S. R. WILLIAMS 
Professor of Physics, Emeritus, Amherst College, Amherst, Massachusetts 


Magnetic Field at a Given Point O Due to a Small Magnet. — Some magneticians are beginning 
to feel that the concept of isolated magnetic poles is obsolete. They insist that magnetic 
problems should be solved in terms of moving electric charges which produce the so-called 
magnetic poles. Until such time as we can express fundamental problems connected with 
magnets in terms of the magnetic moment of a spinning electron, however, we will do well not 
to throw this useful concept of magnetic poles into the discard. 





Coulomb discovered certain laws about the way magnetic poles react on one another and ex- 


pressed them in a single formula with considerable precision. In any medium this familiar 
formula is 


The usefulness of the concept of isolated magnetic poles may be seen in their application to 
the problem of finding numerical values for the magnetic field strengths which surround a 
magnetic dipole. Everyone is familiar with the experiment of scattering iron filings about a 
magnet as shown in Fig. 1. Wherever we see the filings lining up in definite rows, we think 
of them as representing the paths along which the forces due to the poles act. These lines of 
force indicate that we are in the presence of a magnetic field. The magnetic field intensity at 
any point in such a field is, by definition, measured by the force in dynes which that field 
exerts upon a unit magnetic pole placed at the point in the field where we wish to know its 
value. In the following paragraphs we will see how far we can go in solving this problem by 
simple mathematical procedure. 


A magnet, N - S, Fig. 2, has a magnetic field surrounding it something like the one shown in 
Fig. 1. What will be the strength of the magnetic field of this magnet at a point O when r is 
the radius vector connecting O with the center C of the magnet? 6 is the angle which the axis 
of the magnet makes with the radius vector r and f¢ is the distance between the two poles of 
the magnet. In Fig. 2 one may picture the point O as located at the origin of a system of rec- 
tangular co-ordinates and the magnet, N - S, swinging anywhere about this point, but with its 
axis parallel to the Y-axis. When this is done 6 is the angle both between the radius vector r 
and the axis of the dipole and between the radius vector r and the axis of Y. 


At O is placed a unit magnetic pole (positive), because the force on it due to the magnet, N -S, 
will be the field intensity of the magnet at the point O, and the direction in which it tends to 
move will be the direction of the field. This last point is also a matter of definition. 


It will be assumed that this experiment will be carried out in a vacuum, so that the field in- 
tensity at a distance r from a single magnetic pole will be by definition: 


m 
H - 
r 


This follows from Coulomb’s law for the forces between magnetic poles in a vacuum. 
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Fig. 1. The Field of Force 
About a Dipole, 





Fig. 2. The Field Strength of a Magnet, 
N - S, at a Point O, 


In this problem the dipole is treated as though it had two independent isolated magnetic poles. 
It is to be noted that the N-pole can act individually on the unit positive magnetic pole placed 
at O. The same thing holds for the S-pole. The field intensity of the N-pole at O is Hy = m/r* 
and that of the S-pole is Hg = m/r* . 


The total force on the unit magnetic pole at O will then be these two forces indicated above 
and represented in magnitude and direction by the two vectors, On and Os, in Fig. 2. Asa 
matter of convenience it will be found best to resolve these forces into their components 
along the axes of X and Y. When these are found, compound them into their resultant, which 
will be found to be Or= R=Y\Y x’ + Y’, the actual force acting on the unit magnetic pole at O, 
or the field intensity of the dipole, N - S, at the point O, which will be the answer to our 
problem. 


For each pole, N and §S, there will be two components, one along the X-axis and the other 
along the Y-axis, or there will be two components along each axis. Thus, taking the Y-axis 
first there will be the two components along it and the total force acting in that direction will 
be 


sin6d sin8 
F m + .-m—;“ 
 e r 
while for the Y-axis 
—_ cos6, _ - cos8, 








5 r r 
It will now be recalled that in the original statement of the problem the only terms given 
were 8, rand {. In these equations for Fy and Fy there are other values given, all of which 
are unknown. These are 6,, 8,, r,, and r, and equivalents for these must be found in terms of 
the known quantities given. Let us write out all of these relations which we can find and then 
use what we need of them. 


oe 


Digression: r,= r+ re rf{cos@ = r’ - rfcos® 





2 
m-=r f + rf{cos6 = r° + rfcos® 
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when the assumption is made that ( is very small with respect to the value of r. This as- 


sumption will hold throughout this solution. 





r+r,=2r. 1, -7F,= {cos®. r,r, = r’. 
y y- f y+ f 
=cos8. ——— = cos6,. = cos6, 
r r, : 
*=sind. ~-=sind. = = sin6, 
, 1 r, 
{ sin@cosé (sin’é 


sin6, - sin6, cos6, - cos8, = 


r 


Returning to the original equations for F, and Fy, the foregoing values will be substituted 
wherever needed and the equations can then be transformed as follows: 





. mr, sin, mr, sin@, 
x - VS == — - ee 

on 22 

rr, mri ; 
m 2 , . 
— (r} sin6, - r. sin®,) 
4 l 
r 
m {2 a , 2 > , 7 
— (r° + rf{cos8) sin6, — (r° - rf{cos®) sin6, 
r* | J 
m r 2 P : : ! 
—= |r’ (sin®, - sin@,) + rfcos6 (sin®, + sin@.)| 
| 
m r_2 ({sin8 cos8) : 
= in (fsin 2r sin@ cosé | 
J J 
m , . 

= mf [sin8 cos® + 2sin8cos8| 

r 
M . 
= (3 sin® cos6 | 
r 


where M = mf, the magnetic moment of the dipole. Similarly, 
Fy mr, cos®, mr cos@, 


2.2 2.2 
rir, rr, 


m 
= (r, cos@, - r; cos®,) 


r 
m fy/.2 "4 2 ] 
=z | (tm + r{cos8) cos6, - (r° - rf cos@) cos8, 
y * J 
m i 2 , 7 
= |r (cos®, - cos6,) + rf{cos® (cos®, + cos®,) 
| J 
= r’ (- £ sin’*@)+ 2rf cos’6| 
ry * r ‘ 
= [ - sin’® + 2 cos*6) 
r - - 

: s [cos*@ - 1 + 2 cos*6) 
r | J 


5 [3 cos’6 - 1] 
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As indicated above the total force at O due to the dipole will be 


re’ 
w 
Ss 
@®D 
oO 
je) 
n 
®D 
“we 
+ 
| 
a | 
ie) 
n 
we 
®D 
i 
— 
~~» 


= \ (9 sin*® cos’6 + 9 cos*® - 6 cos*@ + 1) 
r 


—-. . 
= |9 cos’@ (sin*@ + cos*6) - 6 cos’*@ + 1 
| } 
fra ~~ 

— (9 cos’® - 6 cos’@ + 1) 

‘r 

M. (3 cos’@ + 1 

r ¥ 


To arrive at this quite simple result we had to make the assumption that { was small with 
respect to r. This is often what we have to do with mathematics in order to get any result at 
all - make some simplifying assumption. If we had not done it in this case, we would have 
had equations beyond our powers to handle with elementary mathematics. 


Since mf = M, the magnetic moment of the dipole, appears in all of these problems, we could 
have started with it in the beginning and thereby have simplified the calculations greatly. 


In the meantime, interpretation of the final results is in order, because mathematics has no 
practical value or significance unless it can be interpreted. When 6 = 90°, cos 90° = O and, 
therefore, the final result becomes 





This holds for the condition that the magnet stands broadside to the point O. O lies on a line 
which is at right angles to the axis of the magnet erected at the center of the magnet. This is 
known as the B tangent position used by Gauss in his method for finding the value of the 
earth’s magnet field. On the other hand, if 6 is O’, cos® = 1 and 


Ft > M \3+1 aM 


In this case, the point O lies on a line which is a continuation of the axis of the magnet. It is 
the A tangent position of Gauss. If the distance r is the same for both cases, the force at O 
in the latter position will be twice that in the former position. 


Mathematics Applied to a Medium Filled with Magnetic Dipoles 








Having found the expression for the magnetic field around a single magnet, or dipole, it is 
quite conceivable that by a process of integration it would be possible to find at any point 
within a medium filled with dipoles - and all turned in the same direction - the resultant 
force of all such dipoles. One may set up such an equation, but its solution is too complex 
for our purposes in dealing with simple mathematics as a tool. 
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By the use of a little imagination, however, the problem can be solved quite readily. The 
following solution also illustrates what has already been mentioned, that engineering and 
physics are rich in pictures, and all that is needed is to find the right picture which can be 
expressed in simple mathematics. 
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Fig. 3. Resultant Magnetic Field Due to 
Paired Dipoles 


In Fig. 3 is shown by a solid line circle the cross-section of a very large sphere (b) of radius 
r, filled with magnetic dipoles turned in helter-skelter fashion. This group of dipoles has a 
magnetic field H applied to it in the direction shown. This turns all of the dipoles into the 
Same direction as the applied field H. In this state the medium is said to be magnetically 
saturated. For simplicity’s sake only four dipoles, Nos. 1, 2, 3, and 4 will be indicated, but 
they will show that when the sphere is magnetized to saturation, the positive or north-seeking 
poles will all lie in the sphere (a) whose cross-section is represented by the circle drawn in 
short dashes, while all the negative or south-seeking poles will be located in the sphere (c) 
represented by the circle drawn in long dashes. 


We now raise the question: What is the force at O, the center of the original sphere, when the 
field H is withdrawn but all dipoles still remain in the position of magnetic saturation and 
each one exerts a field at O given by the equation 


F = \ /3cos?6 + 1 
r ‘ 


If one had a mass m, say near the middle of the earth where the distance to the center was 
one mile, the effective gravitational force acting upon the mass m would be due only to that 
part of the earth which lies within the sphere of the one mile radius. While all the rest of the 
earth would be acting upon the mass m, there would still be a neutralization of all of these 
forces, and only the pull of the sphere of one mile radius would be effective in attracting the 
mass m. A similar situation exists between the unit magnetic pole at the center of the origi- 
nal sphere and the attractions which the positive and negative poles in the other two spheres 
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exert upon it. The unit magnetic pole at the center of the original sphere is off the center of 
the other two spheres by a distance which is equal to half the length of the dipoles. Call this 
distance {/2. Consequently only the positive poles in the upper small sphere (d) will have any 
effect on the unit magnetic pole, m’' = 1, and the same will be true for the negative poles found 
in the small sphere (e) below the center. The result of these forces on the unit pole is that 
the positive poles in the upper sphere push m’ = +1 downward toward the bottom of the page, 
and equally so the negative poles in the lower small sphere attract m’ = +1 downward toward 
the bottom of the page, and we see that the resultant force at O is in a direction opposite to 
that which H originally imposed on the point O. This is why we call this a demagnetizing 
force. The radii of the two small spheres will both be equal to (/2, hence the volume of each 
one of the small spheres will be 


Multiply this volume by the number of magnetic poles per unit volume (= n) to find the number 
of positive, as well as negative, poles in the upper and lower small spheres respectively. The 
total number then will be 


mn 


ae <9 


Now apply Coulomb’s law to the forces between the poles in each of the small spheres and the 
unit pole, m’ = +1. 


-NX1 -nmnf* - 2nnmf 





F, = 





(f)2 6{? 3 
2 4 
eet Nxl -nmnf{* -2nnm/f 
— (fe BP : 
2 4 


Where m, = m, = m, the pole strengths of the dipoles imbedded in the medium. 


The total force at O, therefore, will be 


Fy (-F,) (—F_) = 2F 
-4nnmf -4 
he 3 mnM 


where M is the magnetic moment of each dipole being considered. Inasmuch as n is the num- 
ber of dipoles per unit volume, then nM is the magnetic moment per unit volume of the me- 
dium containing all of these dipoles. By definition, the magnetic moment per unit volume is 
the same as the intensity of magnetization ( and here it is measured within the medium. 
Since nM =I, 


-4nI 
3 





Ft 


This equation gives the value of the force which all the dipoles in the medium contained in the 
large sphere of radius r exert upon a unit magnetic pole at the center of the sphere. The 
results are familiar equations. 
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Critical Angle, 6c = 26. In making the drawing shown in Fig. 3, it was interesting to see how, 





as magnets 1 and 2, moved symmetrically around the circle, i.e., in a plane from 8 = 90° to 
8 = O°. there was a definite line at which, as the centers of the magnets moved over these 
lines at A and C, or E and F, Figs. 4 and 5, a magnetometer needle at O reversed its direc- 
tion of pointing. This was observed by simply placing at O a magnetometer needle and sliding 
the magnets 1 and 2 symmetrically around on the circle but still keeping their constant di- 
rection of pointing. See Figs. 3 and 4. Thus if the two magnets were in the sector between 
lines OA and OC, Figs. 4 and 5, the magnetometer needle would point toward B. The same 
thing happened when the magnets were between the lines OD and OF. On the other hand, if 
the magnets passed over into the sectors between the lines OA and OD, or between OC and 
OF, the magnetometer needle would point toward E. There seemed to be very definite bound- 
ary lines between the four sectors indicated in Figs. 4 and 5. 





109-28 


> 


70 32 


Fig. 4, In the solid cones, ABCO 
and DEFO, the dipoles contained 
therein augment H in the Y di- = 


rection. All dipoles in the sphere 





outside of these two cones op- 

pose H in the Y direction. 
Fig. 5. Ina sphere filled with uniformly dis- 
tributed dipoles, the contribution ofany given 
dipole tothe resultant field at O depends upon 
its position in the sphere. 


In Fig. 5 only one movable dipole is indicated in order that the field surrounding it may the 
more readily be imagined, showing how it operates on the magnetometer at O. In dealing 
with symmetrically placed dipoles the components of the magnetic forces between them and 
parallel to the line joining their centers annul each other. 


From the experience in finding the magnetic field about a single magnet, it seemed quite 
possible to find what the combined fields of two dipoles would be. Thus the study has ad- 
vanced to a consideration of two or paired dipoles. Magnets 1 and 2, Fig. 3, were taken first 
and the question raised: What will be thd field intensity of two such dipoles at a common 
point O, and how does the angle 102 between them affect that result? The preliminary ex- 
periment, outlined above, indicated that there was a certain critical angle, 26 = 2102 8. 
where the effect of the paired dipoles would be equal to zero. 
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Here is a splendid opportunity to show how simple mathematics may be used in telling us a 
great deal about the question which we have just raised. It is only going another step in the 
Same procedure which was used in finding the magnetic field surrounding a single magnet. 
Furthermore, the drawing in Fig. 3 was at first simply by chance, yet as we advanced in 
Studying this problem it didn’t seem sc out of place when one thought of iron atoms being laid 
up in cubic form. Here we are dealing with an octet of dipoles placed at the corners of a 
parallelepiped which is our unit to give zero field at the center of the parallelepiped. A brief 
study of Fig. 3 shows that it is merely a doubling of the problem shown in Fig. 2. Instead of 
dealing with the two poles of a single dipole, here we are handling the four poles of two di- 
poles and the procedure becomes exactly the same as that previously considered. 


Mathematical Determination of the Critical Angle 








Case I. Inasmuch as there will be 4 poles for the pair of magnets which will operate upon the 
unit magnetic pole, m’= +1, at the center of the circle, there will be four components to add 
along the X-axis and four along the Y-axis. 





We start by writing the forces which the poles of the two equal and paired magnets exert on 
m’ = +1. This is done by the use of Coulomb’s law for magnetic poles. For convenience in 
starting, write down in a tabular form the forces due to all four poles participating in the 
study, remembering that we assume that all four poles are equal and may be represented by 
the letter m. Care must be used in writing down these individual forces that the proper signs 
are attached to them so that they will be in harmony with the signs usually assigned to the 
various quadrants used in geometry. 


TABLE 

Force X Component Y Component 
F, = m/r; - msin®, /r- + mcos®, /r* 
iF. =m/r | + msin®, /r? - mcos® ? | 

oo, 2 Sinv,/T, . os | 
'F. = m/r? - msin@,/r? - mcos@,/r 

3 2 a 2 2 2 2 
F, = m/r} + msin@ /r* + mcos®8, /r’ 


4 


It will be observed that when proper signs are given to the force components the total force 
along the X-axis will be F, = O, while along the Y-axis 


F.=(F +8) - (F, + F,) = Cem) - Acom™, 
y 1 4 e 3 r 7 

When this is evaluated we shall have the total resultant force acting upon m’= +1 at O, due to 

the dipoles, 1 and 2. This can be carried out in exactly the same way as the earlier equations 

were carried out. Again we shall have the same unknowns as we did earlier; these need not 

be rewritten here. Putting the unknown values of 6,,6,,r,,and r, in terms of the known values 

of 6, rand ¢, and substituting these quantities in the equation for Fy we get the following form: 


. 2mr.cos8, - 2mricos8, 





re 

Fy “n 72%) (- sin*@) + darcos'8| 
i (- sin’@) + 2cos"@ 
4 | 


= 2cos’8 - sin’@ 
| 


17 











An interpretation of this equation shows that the only condition under which Fy O is when 
the quantity within the squared bracket becomes equal to zero, or 


whence, 


2cos*@ = sin’@ = 1 - cos’@ 
3cos’8 = 1 j 
cos8 ie 0.57735 

@ = 54° - 44’ = 54.733+ degrees 


6. 26 = 109.466+ degrees } 


Case Il. Paired Dipoles 2 and 4, Fig. 3. Of equal interest is it to see what effect the paired 


dipoles 2 and 4 will have on m’ = 1 at O. Again we may write expressions for the forces 
which the four poles exert on m’ = 1 and tabulate the components as we did for the magnets 


1 and 2. 


Again, 


while, 


TABLE 


Forces X-Components | Y-Components | 





2 : 2 
F = m/r, - m cos$,/r{ -m sing,/r, | 
| | 
2 ner 2 . , 2 
F, = m/r, + m cos,/r, m sin@,/r, 
2 att . , 2 
F; = m/r, - m cos§,/T; m sing,/r, 
2 . 2 
F. = m/r, + m cos@,/r} - m sing,/r) 
A i —_ — 
F,=0 
2m sing, 2m sing, 
Fy - 2 a 2 ? 
r2 ry 
2m f 7 


— iri sing, - r’ sing, 
2 | 
r ~ 4 


As in Case I, one must substitute for the values of the unknown factors their equivalents in 
terms of the known terms. See above for these values. In terms of the known, 


— |(r*+ 2rf sin) (sin@,) - (r* - 2rf sinO) ( sing,)| 
r L 


2 ~- _~ 
L. |r*(sing, - sin@,) + 2fr sin@ (sin@, + sin@2) | 


2m 6c :, 9 sa) 
— |r(sin@, - sin@,) + 4{r sin’¢ 


— |2fr(sin’@ - 1) + 4¢r sin’¢ 


<i [sin’@ - 1+ 2sin’@| 
, Ss - 


” (3sin’¢ - 1] 
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Interpreting this equation, it may be seen that Fy O only when the quantity in the squared 
brackets is equal to zero, or 


3sin*@ = 1 
sin’? 
ig 
sing \3 0.57735 
9 = 35° - 16’ = 35.266 degrees 


The results obtained in Case I and in Case II may be thought of as parts of a common picture 
in any plane drawn through the sphere and including the Y-axis. In a simple drawing like Fig. 
4 it may be seen how the angle found in Case I and that in Case II are related to each other. 
It will be seen that they are supplements, one of the other. That these angles are correct, as 
found mathematically, is demonstrated in Fig. 6 for Case I and in Fig. 7 for Case II. 


tht di 


SON by; 
; ; 4 \ Fig. 7. Around a line containing the axes of 
,/ ’ the two magnets is a circle of zero, 
—_— or neutral, points. The plane of this 
circle is normal to the plane of the 
Fig. 6. Neutral points to two similar magnets paper and intersects the latter of the 
when broadside to each other. points A and B, 


Four dipoles laid at the points A, B, C, and D in Fig. 4, or at points 1, 2, 3, and 4 in Fig. 3 
would produce a zero field at O. Furthermore, if we consider ABCDEF as the periphery of 
a large sphere inclosing uniformly distributed dipoles throughout its volume, it may be seen 
that all dipoles lying in the two solid cones, ABCO and DEFO, within the sphere, give a re- 
sultant field at O which is directed in the same direction as the applied field H, while all di- 
poles lying in the band which encircles the sphere, represented in cross-section by the tri- 
angles, AOD and COF, give a resultant field at O which is opposed to H. Here we see why the 
conical poles of an electromagnet are the most efficient when ground to an angle of 109.46°. 


It is beyond the limits of this paper to go into a discussion of how a magnet shaped like the 
one shown in Fig. 7 would behave. Would it have a greater magnetic moment per unit volume 
than any other shaped magnet? 


By simple processes of integration, it may be shown that the volume of the two solid cones, 


ABCO and DEFO form 0.42262 part of the total sphere and that the belt around the sphere, 
represented by AOD and COF, form the remaining, or the 0.57738 part of the sphere. This 
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Fig. 8. Cross-sectional View of 
“El Diablo” magnet Fig. 9. Ahollow plastic sphere through 


which strings have been pulled 
to outline the solid cones, 


ABCO and DEFO, 


means that there are more dipoles helping to oppose the field H at O than there are co-oper- 
ating with it. Fig. 9 is an attempt to show in 3-dimensional space the two solid cones, ABCO 
and DEFO, and the belt around the sphere represented in cross-section by the two triangles, 
AOD and COF. For lack of space the physical interpretation of the mathematical results 
must be reserved for another paper in some physical periodical. 


By some very simple mathematics we have arrived at definite results which could be verified 
experimentally, and have thereby demonstrated that simple mathematics can be a valuable 
tool in science. 
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HIGH-EFFICIENCY CANTILEVER SPRINGS 


DOUGLAS G. ANDERSON 
McIntosh Stamping Corporation 
Detroit, Michigan 


Introduction 


Springs made by a stamping process from flat sheet material are usually rectangular in 
cross-section, although they are frequently bent or curved in profile. When material cost is 
high and the volume of production large, it would often be economical to make the increased 
investment required for tools to make a more efficient form of spring. So far as the author 
knows, however, there has been no practical method for selecting and specifying such springs. 


Nomenclature 


The following nomenclature is used: 


A Area of cross-section normal to stress. 

c Maximum value of y for a given cross-section. 

© Ratio of width to thickness (for Part II of the section shown in Figure 1). 

E Modulus of elasticity. 

e Efficiency. 

eg Gross efficiency. 

f(k) Function of k appearing in the expression for Z. 

F(k) Function of k appearing in the expression for Q. 

I Moment of inertia of cross-section. 

k Ratio of depth of forming to thickness (for section shown in Figure 1). 

L Length as measured along neutral surface from point of application of load. 

M Moment. 

Pp Load. 

Q A certain function of the proportions of a spring. 

R Radius of curvature of neutral surface. 

S Stress. (When not under an integral sign, S is the maximum stress acting in the region 
considered). 

Sa Highest stress allowable for the application. 

t Thickness of sheet material. 

U Energy of deformation of an elastic body. 

V Volume. 


Vg. Gross volume of untrimmed blank. 

W Width of section. 

Wp Width of blank required to produce a section of width W and depth (k + 1)t. 
Ww, Ratio of (Wp - W) tot, a function of k. 


y Distance from neutral axis of cross-section. 

Z Section modulus, I/c. 

fa Deflection at the point of application of load. 

5 Deflection of edge of section with respect to neutral axis, due to distortion of the section 


(see Figure 5). 


Efficiency 


The essential function of a spring is to deflect under load and so to store energy. The energy 
stored in an elastic body of uniform cross-section, when it is uniformly and fully stressed, is 


1 1 1 S._°V 
= e = “a > 2 e See 
U 5 PO 5 (SgA) (+ Ll) +33 AL = 4. 
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If the stress is not uniform throughout the body, the energy stored will be less. The strength 


of the material is then not fully utilized. With respect to energy storage, the efficiency of the 
body is 


which holds for a body of any shape. 


For economical tooling and handling costs, usually a spring will be blanked from a strip just 
wide enough for one trimmed blank. The volume of a length of this strip just long enough for 
one trimmed blank is then the amount of material actually used to make the spring (neglecting 
the value of the scrap). This volume is called the “gross volume,” and the gross efficiency is 
defined as 


Choice of Section 


The form of cross-section to be used should be symmetrical about an axis at right angles to 
the neutral axis to avoid deflection normal to the direction of application of the load. It should 
concentrate as much of the metal as possible as far from the neutral axis as possible to 
maximize J,S dA. It should be readily formable from material of uniform thickness. These 
considerations lead to the form shown in Figure 1. Its properties, and in particular its sec- 
tion modulus, are functions of t,k, and C. The theoretically preferable sharp corners and 
vertical walls are modified for ease of forming in dies without close clearances. 
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Load and Deflection 


The problem is to choose t and to vary k or C or both along the length of the spring in such 
a way as to obtain high efficiency with the load and deflection required. The first step is to 
obtain a relation involving the deflection. 


Neglecting shear stress, the energy equation for a beam with a single concentrated load is 
i l 
2U = PA = f-—— E fu fa S’dAdL (1) 


Assuming, as is usual in elementary beam theory, that the stress at any point of a section 
taken normal to the neutral surface is proportional to its distance from the neutral axis of 
the section, we may write 
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Pa a5. m [SaydA] aL (2) 


where S is the greatest stress within the area A. But the integral within brackets is, by defi- 
nition, the moment of inertia of the section. Substituting in Eq (2) and dividing by P, we have 
the deflection formula 


SI 
4 =—— S, =a (3) 


Eq (3) can easily be derived from the deflection formulae of the curvature method, Van den 
Broek’s elastic energy method, or the theorem of Castigliano. It is unusual in that I appears 
in the numerator of the integrand. This simplifies the evaluation of the integral in many 
cases where I is a function of L. 


At any cross-section of a cantilever we have 


M= PL s < SZ (4) 


For the highest efficiency, S should obviously be independent of L, that is, the spring should 


be a beam of constant strength. Then dL 5 dZ and Eq (3) may be written 
s° 7, 


for a cantilever beam of constant strength. 
Returning to the consideration of the section shown in Figure 1, we now wish to find its section 
modulus. It is convenient to consider parts I and II separately. By definition, I = I, = I,,. But 


also, C; = Cy; = c, and therefore 


I 
: iu Z,+ Zy (6) 
Expression Z, and Z,, as functions t, k, andC, we have 


Z = t* [fy(k) +C fy(k)] (7) 


The function fj is not readily expressible in terms of k directly, but Table I gives values for 
both f, and f;; for values of k from 0 to 10. From Eqs (4) and (7), we can write 


sar = fy(k) + C fyy(lo (8) 


This equation has been plotted for several values of C in Figure 2. 

It is impractical to extend the constant-stress condition to the region where k is less than 1. 
A smooth variation in depth of section is obtained by making k proportional to L from k = 0 to 
k = 1. 

Although in some cases it might be advantageous to vary C with L, particularly in the region 


where k is small, we will confine our attention to springs of uniform width, where C is con- 
stant. Then each curve of Figure 2 represents a member of our family of springs. 
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Substituting Eq (7) directly in Eq (5) and expanding we have, since c = lik + 1)t, 


2 
s*t° 2 
6= Be Se wee [fh + Chyle)- 4 [1,00 + £00] 
_2s't® 


{yr 4,0) dfy() , [s, fy(k)dfyy(k) peat c2 (, Lyp(k) dt (b) 
PE | k¥1 k+1 k+l k+1 


Since the integrals are functions of k alone, we can write 


are Fo(k) + C F,(k) + C* F,(k) = Q(c,k) (9) 


Note that Q is a quadratic in C. Our interest in it stems from the condition imposed upon it 
by the deflection. However, it has a more general significance than the method of derivation 
might lead one to suspect. It might be called a “shape factor”; by the choice of suitable 
parameters in place of C and k it is definable for any class of homogeneous elastic bodies. 


For a cantilever of uniform rectangular section, Na aq). For a constant-strength 
cantilever of constant thickness, arr age). 


Values of F,, F,, and F, are given in Table I. The values for k = 1 are the coefficients of a 
quadratic in C which, for values of C from 4 to 32, closely approximates the exact solution 
for a beam in which k is proportional to L from k= 0 to k=1. The values of the F’s for 
higher values of k were obtained by numerical integration. 


Table I also gives values of w,, a function of k needed for determining the width of the flat 
blank from which the spring is to be formed. At any section, the width of blank is 



































Wp = W + wt (10) 
TABLE I 
[xe [oyfTiy2]3 [4 [5 | 6 7 Te | 9 | 101 
as a —+———— - — a —_—_—_—+ — + $$$} 

i(k) |. 33 | 1.56 |3.29 | 6. 02 | 9.55 | 13.83 | 18.84 | 24.55 | 30.75 |37.70 | 45.8 
- + ——— ee — Tr 
f;(k) | .167 | 333 | .722 [h. 167 i 634 | 2.091 | 2.595 | 3.083 | 3.574 | 4.067 | 4.560 
-—— + + - + - +-——_—- Tana, ana loam & fo 

F,(k) | 0 |.75 [4.2 | 11.2 | 23.4] 41.5 | 66.6 | 99.6 [139.9 | 1189.9 |246.3 
He + + + : —_— — 

F(k) | 0 |.58 |2.06 | 4.66 | 8.45 | 13. 38 | 19. 43 | 26.56 [34.60 | 43.73 [53.48 
'F,(k) | 0 |.054|.212 | .449 | .737 | 1.061. T1411 11.780 [2. 164 Tas 560 | 2.965 
+ + + + —* + t — — —+— 
fw. | o | 1 | 7 [18 | 33 | 5.0 | 6.8 | 8.7 | 10.7, [12 12.7 | 14.7] 





Eq (10) was plotted against k, for values of C from 4 to 32, and the lines of constant Q were 
transferred to Figure 2. 


Use of the Chart 





Suppose that the spring is required to have the following characteristics: 


P = 100 lbs. QO = .25 in. S = 100,000 lbs./sq. in. 
E = 30,000,000 lbs./sq. in. W= .96 in. 





Assuming t 


.030, solve for Q. 














OPE .25x10*x3x 10’ 
Q* SF * ox x 107 ~ 9085. 
From Figure 1, we have 
Cc x - 8 (11) 
Solving for C Cc -+ 2 32 - 8 = 24 
tin dint. 030 ~ : 
At the point (Q = 3085, C = 24) of Figure 2, we read that k = 9.6 approx., and ey 148. 
43 - 5 3 -6 } 
Solving for L, L a =; <n eS | aa 4.02. 
From Figure 3 we read, at k = 9.6, wy, = 13.9. 
Solving for Wp, Wp= Wt wyet= .96+ 13.9 x .03 = 1.38. 


Intermediate values of k and Wp are easily determined. 


= . With this value, and C 


For any desired value of L, compute 


24, find k from Figure 2, then w, from Figure 3, and compute 


Wy. These data, arranged in tabular form, will be used by the draftsman and the diemaker. 


Efficiency and Distortion 





In the above example, t was chosen arbitrarily. 


A larger value would yield smaller values 


for C and k. Returning to our expression for gross efficiency, and substituting in terms of 


our new functions, we have 
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Figure 4 


Figure 4, which is a graph of Eq (12), shows that there is an optimum value of k for each 


value of C. but it also indicates that large values of C are desirable. 


This judgment must be 


modified by consideration of distortion and possible buckling. 


Curvature of the spring due to its own deflection causes a distortion of each cross-section, 
of the type shown in Figure 5. The value of 6 is given approximately by 
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eh es hs 3G)" + ye? + 1 [685 + 64(9) + 36(9)° + 24(6)° | } (13) 


This formula is based on the assumption that the neutral surface is initially plane, and that 
§/t is so small compared to k that the resulting change in stress distribution may be neg- 
lected. It does not apply when k is less than 1. In general, this distortion is most serious 
whenk 1. Figure 6 shows 6/t as a function of C when k= 1 and S/E = .0033, as in our ex- 
ample above. Obviously it would be inadvisable to use a value of C much greater than 24. 
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Figure 5 Figure 6 


If the neutral surface of the spring has appreciable curvature in the stressed condition, the 
limitation on C is more severe. Let R be the radius of curvature of the profile view of the 
neutral surface in the stressed condition, positive when the center of curvature is on the 
compression side of the spring. Then the factors outside the curved brackets on the right- 
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hand side of Eq (13) are replaced by a . 


It is favorable to stability for the curvature to be negative, so that I increases with the dis- 
tortion, and for the free edges of the spring to be on the tension side, to reduce any tendency 


to buckle by column action. 


Comparison with Conventional Springs 





Where design conditions allow the use of this ribbed spring, its advantages are truly remark- 
able. Table II compares the spring of our example with the triangular and uniform-section 
flat springs which meet the same design requirements. The material saved is from 50 to 75 
percent. The stability with respect to variation in stock thickness is often a prime considera- 
tion. An advantage not shown by the table is the greatly increased accuracy of profile made 
possible by the large depth of section, compared to stock thickness. The author has also had 
occasion to use aluminum for a spring of this type, in an application where endurance limit 
was not a factor, but where resistance to weathering was very important. Replacement of 
the steel spring previously used, which was stressed to 100,000 psi, by an aluminum spring 
stressed to 40,000 psi working in the same space, was possible only because of the additional 
degree of freedom allowed by ribbing the spring. The savings due to elimination of heat- 
treatment and plating costs, and the higher efficiency of the shape, made the new spring less 
expensive than the one it replaced. 


TABLE I 

= p Aaa Gross | 2changeinrate  _—_ 
Type of Spring | Length | volume | “change in stock thickness | 

| | [Stress |" Deflection | 

Ribbed | 4.02 | 166 | 1+ | 1+ | 
_ - +- +— - —— - + ————_4 
Uniform Rectangular Section | 4.29 | 337 2 3 
re —-4. + -- t ~ 4 
| Triangular Flat | 3.06 | .675 2 | 3 
L " i - = - 1 4 . oumenil 
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AN APPLICATION OF LIEBMAN’S ITERATION METHOD TO VISCOUS FLOW 


LEONARD S. LUSTICK 
Research Laboratories Division, General Motors Corporation 
Detroit, Michigan 


Introduction 





The solution of many problems arising in applied mathematics entails solving a particular 
partial differential equation with its associated boundary conditions. The nature of the differ- 
ential equation and its boundary conditions makes it often impossible to attain an analytical 
solution to the problem. The Liebman Iteration Method is a numerical method of attaining the 
solution at a discreet number of points. 


It might be well to reflect a moment on the sifnificance of a differential equation and how 
such an impersonal sign language as mathematics can represent and predict the behavior of 
an actual physical situation. 


The differential equation itself is only a mathematical statement which is consistent with 
basic physical laws, such as the familiar F = Ma, the conservation of matter and energy, etc. 
Certain assumptions are made in writing the differential equation. How well the equation 
represents the physical situation is dependent on the accuracy of the assumptions, which is in 
turn dependent on the observer’s knowledge or feeling of the physical situation 


The partial differential equation by itself has an infinite number of solutions. However, the 
solution we are interested in is determined uniquely by the singular boundary conditions. 
This last statement is not obvious, and would require a considerable amount of mathematical 
maturity to prove. Since this maturity is not at the disposal of this writer, a bit of mathe- 
matical philosophy will be substituted for it. This philosophy is known as Hadamar’s prin- 
ciple and, in effect, states that if a governing partial differential equation arises out of a 
reasonable physical situation, and it is satisfied everywhere in a particular domain, including 
the boundaries, the solution thus obtained is unique. 


Statement of Problem 





In the problem concerned with in this paper, it is desired to determine the outlet pressure of 
a pump necessary to support a load W on a stationary cylindrical shaft and to determine also 
the pressure distribution under the shaft. The geometry of the problem is shown schemati- 
cally in Figure 1. 


Since the solution obtained by Liebman’s method applies only to a specific geometry, it will 
be necessary to assign dimensions to the pertinent parameters of the problem. These param- 
eters will have the following values: 


Diameter of inlet hole, d = 1/4’’ 
Length of shaft, 1= 1-3/4'’ 
Diameter of shaft, D = 1-1/2’’ 


Applications for this type of problem might arise where it is necessary to avoid the initial 
boundary lubrication conditions upon starting and before the shaft has accelerated to speed. 


28 











It might also be used in problems where the relative motion between the shaft and bearing is 
negligible and no hydrodynamic oil film can be sustained, such as in the wrist pin of an auto- 
mobile. 
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Basic Equations and Assumptions Used in Deriving the Governing Differential Equation 





The basic physical laws used in deriving the governing differential equation for the above 
problem are the familiar F = MA and the conservation of mass. The utilization of these laws 
; , , 1 , — CP. FP. neat 

in the light of rational assumptions results in the governing equation, Be S$ 0, which is 


recognized as LaPlace’s equation with pressure as the dependent variable. 
The assumptions used in obtaining this equation are listed below: 

) The acceleration of the oil is zero. 

) The effect of curvature is neglected. 
) 
) 


Steady-state conditions exist. 
The oil is an incompressible fluid. 








Having once determined the governing equation, it is advisable for purposes of the iteration 

method to reduce it to a non-dimensional form. In our particular case this can be done quite 

_ O&P, &P 

easily by multiplying the equation, $x * Sy 
. ' 

: P 

a + tS o, where U equals the non-dimensional parameter P, 


o by 5. The equation is now reduced to 
oO 








Reduction of Governing Equation to a Difference Equation 





With the equation in non-dimensional form, the next step is to reduce it to a difference equa- 
tion by means of finite difference theory. This difference equation can be realized in the 
following way. Suppose for purposes of illustration we visualize the domain of our problem 
divided into a square mesh as shown in Figure. 2. U, is the value of the function at one arbi- 
trarily chosen discreet point, O, in the mesh. 


U., U,, U,, U, are values of the function at the four surrounding points. Then it follows that 
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In like manner, by —— a? 
Hence, the governing differential equation o o'u 0, can be written for arbitrary point 0 
— . a P a 
; g £ q ‘Sx by y Pp 
as the difference equation, U. + U,+ U, + U, - 4w,= 0. 


At points adjacent to the boundary it often happens that the surrounding points are not equally 
spaced around it, as shown in Figure 3. In the event of this condition, the difference equation 
at this type point is 
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The problem of iteration is at once obvious and is seen to consist of simultaneously satisfying 
the difference equation at a large number of specified discreet points in the domain. 


The final result of the iteration procedure is shown in Figure 3. 


The value of P, (the inlet pump pressure) was obtained by numerical integration, taking cog- 
nizance of the fact that the bearing is cylindrical and adding only the vertical components of 
the pressure. These, of course, must be equal to W. 


Summary of Iteration Method 





(1) 
(2) 
(3) 
(4) 


(5) 


(6) q 


Derive from basic physical laws the governing partial differential equation. 

Non-dimensionalize the above equation. 

Transform the equation in (2) into a difference equation. 

Establish the mesh and make original assumption as to the value of the function at the 

nodal points. This original assumption can save considerable time if accurate. 

Proceed by an orderly method to simultaneously satisfy the difference equation at the 

nodal points. 

Check to see if an advance to a finer net is required. Referring to Figure 2, the pro- 

cedure is as follows: 

(a) W is calculated as average of A, B, D and E. X, Y and Z are calculated in like 
manner. 

(b) E is then calculated as the average of W, X, Y and Z. 

(c) The above value of E should check closely with the previously determined value. If 
this condition is not met, it is an indication to advance to a finer net using the pre- 
vious approximation and linear interpolation to determine the original assumption 
the finer net. 


In closing, it should be pointed out that no attempt at originality has been made aside from the 
specific geometry of the problem which is similar to a problem solved in a graduate course 
at Syracuse University. 


The purpose of this paper has been to illustrate the steps in using the Iteration Method in 
solving a boundary value problem. 
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PERFORMANCE OF CENTRIFUGAL PUMPS FROM AFFINITY LAWS 
AT NON-CORRESPONDING OPERATING POINTS 





ARVID E. ROACH and HOWARD H. KEHRL 
Research Laboratories Division, General Motors Corporation 


Detroit, Michigan 


Introduction 


This paper is concerned with the problem of governing a variable-speed centrifugal pump 
drive. The problem is: In what manner shall the speed of the pump be varied in order that 
its head-capacity curve may intersect the head-capacity curve of the receiver system at any 
specified point on the latter curve? The solution proposed in the paper is illustrative of a 
method of handling such problems by means of the principle of dynamic similitude, which 
method may be useful where not enough is known about the physical situation to enable one to 
solve the differential equations for the system. 


Nomenclature 


The following nomenclature is used: 


Symbol Dimensions 
~ Q Capacity, or quantity of fluid passing through the pump, cubic feet Pe 
per second 
H Head developed by pump, or energy transferred from pump to fluid, Lt 
foot-pounds per slug 
N Speed of pump impeller, revolutions per second a 
D Diameter of pump impeller or other characteristic linear dimension,feet L 
p Density of fluid flowing, slugs per cubic foot ML~°* 
ub Viscosity of fluid flowing, slugs per foot-second ML~T™ 
r Ratio of any linear dimensions of pump Dimensionless 


Subscripts 


1 Initial conditions 
2 Final conditions 

m Model pump 

p Prototype pump 


The plot of the variation in head developed by a centrifugal pump with variation in capacity at 
constant impeller speed is called the characteristic curve (Fig.1) of the pump. This curve is 
obtained experimentally, usually at a single impeller speed. The variation of head and ca- 
pacity at other speeds retains the general form of Fig. 1 and may readily be computed by 
means of the affinity laws for centrifugal pumps. 


The affinity laws were first found experimentally but may be derived by dimensional analysis. 
From experimental and theoretical considerations, it may be shown that the capacity cf a 
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pump is a function of the conditions of operation, the size and proportions of the pump, and 
the properties of the pumped fluid, that is, 


Q= f{(H, N, D,p,p,r’,r’’,.. . r®) (1) 


where r’,r’,...r® 


are such ratios of linear dimensions of the pump as may be needed to 
describe its geometry. 
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By the nm theorem, Eq (1) may be rewritten in the form 
m, = F(t, M5, . . . %-5) (2) 


where 7,,1,,...,_, are the n - 3 possible independent dimensionless ratios which may be 
formed from the variables of Eq (1). 


Sclving for the dimensionless m’s, we obtain 


3 


m,= QN” D~’*, 1,= HN’ D”, 1,=pND up , n, =r’, etc. 


Eq (2) becomes 


Q=ND°- FKHN” D”’,pND p”,r',r’,... 3") (2a) 


where the form of the function F is unknown. 


Let the subscripts p and m designate a prototype pump and a model pump respectively. Then 
Ny D,’ - F(Hp Np~ Dp’, PNp Dp pf ', fp’ To's --- > Fp 

Sp = F r, P P Zz wy p Ms nf Pp p Th (3) 
Qm Nm Dm: Fim Nm™~ Dm > PpNmDm } > lm’, Tm’;--+>lm ? 


Eq (3) provides a basis for comparing pumps, that is, for predicting the performance of a 
prototype pump from the known performance of a model pump whenever the separate argu- 
ments of the unknown function F are equal, that is, whenever (1,)), = (T)m, (%s)p = (M3)m, ete. 
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The parameter 1, = pN D’ uw is Reynolds’ number. Unless different fluids are used in the 
prototype and the model, it will in general be impossible to make (N,)p = (,)m, while keeping 
(M2)p (T1,)m, etc. However, 1, affects primarily the fluid friction of the pump and, as shown 
by Fig. 2, large changes in Reynolds’ number have relatively little influence on the hydraulic 
resistance to flow so long as the flow is fully turbulent. The influence of Reynolds’ number 
on the operation of centrifugal pumps may therefore be assumed negligible, an assumption 
which is justified in practice by the virtual independence of the hydraulic efficiency of centri- 
fugal pumps upon Reynolds’ number. 





Figure 2 


The conditions (m,) (T)m, etc., may be fulfilled simply by making the model a geometric 
image of the prototype, a condition which is easily satisfied on a macroscopic scale. 


The affinity laws thus reduce to the statement that 


Np Dp’ H Ny,’ D,’ 
Sp _ when P P,P 


- (4) 
Qm Nm Dm Hy Nm Dm”’ 


and vice versa. These laws may be used to investigate the effect of changes in scale of geo- 
metrically similar pumps or of changes in speed for a given pump. 


A simple problem will serve to illustrate some of the properties of these laws. Let it be re- 
quired to determine the characteristic curves for the pump of Fig. 1 at several other speeds 
of operation. Take as the basis for computation the point A having the coordinates Ha = 44, 
Qa = 200, and Na = 1000. At speed Np = 1200, we have 


Qp -Q, NB PB ® - 200 x 1.2 = 240 


Na DA 
r 2 2 
and Hp - Ha xe an 44 x 1.44 = 63.4 


We may tabulate these results with others computed in this manner for different speeds: 


N H Q 
800 28.2 160 
900 35.6 180 

1000 44.0 200 

1100 53.4 220 

1200 63.4 240 





Points computed in this manner are called corresponding points and are characterized by the 
property that n, = constant, n, = constant, etc. Similar groups of corresponding points may 
be computed for other operating points from the given curve, and in this way a family of char- 
acteristic curves may be plotted, as in Fig. 3. Note that all corresponding points lie on 
Square parabolas, such as O-A-B, a fact which is evident from inspection of Eq (4) when 
Dp = Dm. Curves computed in this manner have been shown to correspond very closely to 
experimentally determined characteristic curves. The principal restriction is that one is 
compelled always to deal with corresponding points. 
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Performance at Non-Corresponding Operating Points 








Problems arise in which it is desired to compute operating data which does not lie on the 
parabolic locus of corresponding points. For example, consider a pump delivering fluid to a 
discharge system composed of an arrangement of distribution pipes. The characteristics of 
the pump and discharge system may be such as is represented in Fig. 4. The pump will op- 
erate at point A where its characteristic curve intersects the system curve. Suppose. that the 
pump is driven by means of a variable-speed motor and that it is desired to determine the 
speed at which it must run in order that its characteristic curve will intersect the system 
curve at some other point C not lying on the locus of corresponding points passing through 
point A. The affinity laws, as stated, offer no means of proceeding from point A to point C, 
for these points are non-corresponding. In practice, such a problem is usually solved graph- 
ically by sketching a family of pump characteristic curves at different speeds and by picking 
off the speed of the characteristic curve passing through point C. 


An ingenious short cut, suggested by Mr. John G. Wright (Mechanical Engineering Depart- 
ment, College of Engineering, Wayne University), makes use of the fact that at some point on 
the characteristic curve there is a point D which is in correspondence with the point C. This 
point may be determined graphically by plotting the parabolic locus of corresponding points 
passing through point C. Having determined the coordinates of point D in this manner, the 
remaining coordinates of point C may be found from the affinity laws. 


A graphical solution is of limited usefulness in specifying what kind of governor is needed to 
regulate the pump speed so as to satisfy changing system requirements. What is needed for 
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this purpose is an analytical relation between the speed change 4 N required for a given ca- 
pacity change 4 Q and head change 4 H. 
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Figure 4 


The complete characteristic of a centrifugal pump may be represented by a surface in H-Q-N 
space, as shown in Fig. 5. The system curve is independent of pump speed and may likewise 
be represented as a surface in this space. The operating line of the combined pump and dis- 
charge system will be determined by the intersection of these surfaces. 


The known data of the problem are, for the pump, H,, Q,, N, and xe . These data are not 
L ji 


sufficient to characterize completely the characteristic surface. However, if some other 
fo Q) ; ; ' 

derivative, say a , were known sufficient data would be available to characterize a plane 
‘ | 


tangent to the characteristic surface at the point H,, Q,,N,. The intersection of this plane 


with the system surface would represent, for small changes AQ and AH, the operating line 
of the pump-system combination. 


The derivative ae may be obtained by differentiating Eq (2a) which may be rewritten for 
s j i 


convenience in the form 
QN™D~ = f (HN7~*D™, and other n’s) (5) 


where all the m’s are considered to be constant, except QN~*D~* and HN~’D~ which are per- 
mitted to undergo small changes. The difficulty which arises in differentiating a function 
whose form is unknown has been ably treated by Mr. M. D. Hersey in a paper (“A Relation 
Connecting the Derivatives of Physical Quantities”) published in the Scientific Papers of the 
Bureau of Standards, Vol. 15, 1919-20. This paper, which depends upon and is a corollary to 
Buckingham’s n-theorem, has apparently been neglected by later investigators in this field 
and, in our opinion, merits further consideration. 
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Following Hersey’s method in differentiating Eq (5), we have 























3(QN™*D-*) df _ 0(HN™“*D~*) 
ON d(HN~D~) oN 
df _ -2(HN~*D™~) 6) 
d(HN~D~) N 
But also 
d(QN""D™) _ dQ (QN™"D™) _ (QN=*D~*) (7) 
ON ON Q N : 
From Eqs (6) and (7), we have 
(HN~*D~*) df =» _— -N(QN™'D™*) /1 0Q | 1) (8) 
d(HN~’D~*) 2 \Q ON N/ 
Similarly, we may write 
(HN™*D™) ___df__ = H(QN™'D™) (1 0Q) (9) 
d(HN~D~ \Q in 
From Eqs (8) and (9), we have 
0Q_Q ,H 8Q 
INN ?N OO " 


which is the desired derivative. 
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Since Q = f(H, N), then 


dQ 02 4 + 28 an 





OH ON 
and dN = dQ - ie dH 
—— 
ON 
or dN = dQ - 23 dH 
Q fH 5Q (11) 
N N 5H 


To illustrate the application of Eq (11), we may consider the following problem: A pump sup- 
plies fluid to a system having the combined characteristics shown in Fig. 4. It is required to 
change the pump speed so as to supply a larger quantity of fluid to the system. At the initial 
condition, the pump and system curves intersect at Q = 212.4 gpm, H = 43.2 ft, and N = 1000 
rpm; at this point the measured value of 5Q/5H = 13.65. At the final condition, the system 
curve and new pump curve are to intersect at Q = 220 gpm, H = 45.0 ft. The new speed for 
the pump will be 


7.6 + (13.65 x 1.8) 
N = 1000 + 193 5x43. x13 65 1000 + 23.2 = 1023.2 rpm. 


412.8 | 
0° T1000 





100 


A graphical solution using Wright’s method gives Q = 215 gpm and H = 43 ft for the point on 
the original pump curve corresponding to the final point Q = 220 gpm and H = 45 ft. From the 
affinity laws, N = 1023 rpm. 


Such excellent agreement will be obtained only when the change in operating point is fairly 
small. As the change becomes larger, Eq (11) will yield results which are increasingly in 
error. 
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ANALYTICAL DERIVATION OF GEAR TOOTH ‘‘X’’ FACTOR 


DOUGLAS HUGHSON 
Modern Engineering Service 
Berkley, Michigan 


The tooth form factor plays an important part in most gear tooth bending stress formulae in 
use today. This form factor was originated by Wilfred Lewis and presented in a paper read 
before the Engineers Club of Philadelphia, on October 15, 1892. For standard gears, the 
factor y which is equal to 2X/3p has been tabulated in many machine design texts. These 
tables are of little use in the design of special gears which have become widespread today. 


For years it has been common practice in automotive gear work to lay out the gear tooth ten 
or twenty times size and to scale the X Factor. The accuracy of this time-consuming method 
is of course limited by the accuracy of the draftsman. 


The analytical method discussed here, although theoretically correct only for gears hobbed 
with a sharp corner tool, may be applied even when the hob corners are rounded slightly. 
The factor becomes more conservative as the radius of the hob tooth corner increases. 


For symbols and equations, refer to Figure 1. 


When the tooth center line coincides with the y axis, and the gear center falls on the origin, 
the equation of the fillet becomes 


x =Rsina-bsina-Racosa+ Kcosa (1) 


y= Rcosa-bcosa+Rasina-Ksina (2) 


where R is the pitch radius of gear, K is 5 +b tan By, b is the tooth dedendum, t is the tooth 


thickness at R, and o, is the pressure angle at R. 
Selecting any point x, y, on the tooth fillet, we have 


x 


Ky =A (3) 


where A is the greatest distance from the center of the gear to a point on the tooth center 
line through which the tooth force acts. 


The tooth X factor is defined at the minimum value of Xy: 





dx gy 
dXy _ 2(A-y) x da+ x da _ O 
da (A-y)* 
dx dy 
, qa? £2 4 
or 2(A-y) ZO Xa °° (4) 
dx ' ; 
From (1), = R asine a- K sine a- bcos a= y-Rcosa (5) 
From (2), ay Racos a- Kcos a+ b sinea= -x+Rsinea (6) 
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Then substituting (5) and (6) in (4), we have 


2(A-y) (y-R cosa) - x(x-R sina) = O 
or 2Ay-2AR cosa - 2y’+ 2Ry cosa - x’ + xR sina= O 


When Xy is a minimum, referring to the point on the fillet as x, y, and the corresponding point 
at which the hob pitch line is tangent to the gear pitch circle as x, y,, we have 


2Ay, - 2Ay, - 2y, *+2y, y. - % “+x, x= O 


or 2(A-y,) (y, - y2.) - X, (4, - X%) =M=O (7) 
where y, = R;, cosa + (Ra - K) sind (8) 
X R,. sing - (Ra - K) cosa (9) 
y, = R cosa (10) 
x, = R sineq (11) 
- = (12) 

cos - Bo) 


and rp is the base radius of gear, ®, is the pressure angle at outside radius, B, is the angle 
subtended at center of gear by one-half outside tooth thickness, and R, is the root radius of 
gear. 


Equation (7) may be solved by substituting various trial values of a in (8), and (9), (10) and 
(11) and by substituting values of x,, y,, x, and y, obtained together with ‘A’ back into (7) until 
a value of qa which gives a solution for (7) is obtained. 


By substituting the values of x, and y, corresponding to this value of a in (3), the Gear Tooth 
X Factor is given by: 


x ? 


x : (13) 
A-y, 








Figure | 
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From the foregoing we see that the solution in part consists of a set of equations, namely (7), 
(8), (9), (10), (11) and (12), that must be solved by trial. This theoretically exact although 
cumbersome solution is presented in the hope that someone may discover a suitable substitu- 


tion to simplify these equations and arrive at an analytical approximate solution of sufficient 
accuracy for stress analysis. 


In Figure | a badly undercut tooth has been purposely chosen to show that the formulae apply 
to hobbed involute teeth whether undercut or not 


A sample calculation is detailed below to show the steps to be followed in the solution of the 
problem. 


Mention should be made that Mr. Leonard G. Johnson of the Research Laboratories Division 
of General Motors Corporation has arrived at similar results from a somewhat different ap- 
proach. His results were given in a talk before the Michigan Section of the Mathematical As- 
sociation of America on March 25, 1950, at Ann Arbor, Michigan. 


Sample Calculation 











N= 17 ty = 20° ty .084260 
p-8 2R = 2.125000 A = 1.161698 
b = .144625 2rp = 1.996847 K = .150814 
t = .196350 2Ro = 2.375000 Ry = .917875 
0 ] Xx, y, Xx, Y2 M | 
20° + + + + — , - | 
3400686 107135 937788 363396 998423 000301 | 
19.90° 107242 937343 ( .3161653 999056 000408 
347321 ' ; 
0 
19.9576 107180 937599 362657 998692 000000 
348326 
- L i = | 2 ae 
.1071807 
X = S57999 * -051261 
2x .051261 _ 
3 x .392699 = .087023 
Addendum 


Subsequent to the presentation of the paper, it was brought to the author’s attention that, in 
addition to the paper by Mr. Leonard G. Johnson referred to above, there is a third paper on 
this subject. This is Mr. M. Maletz’ paper on “Analytical Determination of the Form Factor 
in the Beam Formula for a Tooth,” presented before the American Gear Manufacturers As- 
sociation in May, 1941. The papers by Mr. Johnson and by the author deal with the special 
case of a rack-shaped and sharply cornered tool; Mr. Maletz’ paper deals with disk-type as 
well as rack-type tools, having sharp or round or chamfered tips. 
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SURFACE STRESS AND STRAIN OF SPIRAL INVOLUTE GEARS 


VICTOR E, FRANCIS 
National Broach and Machine Company 


Detroit, Michigan 
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The design of spiral or crossed-axes gears has always been a difficult problem for several 
reasons. 


The first part of the design problem is to obtain the proper helix angles of the gears. There 
are an infinite number of usable helix angles for any set, given the center distance and the 
desired ratio. Each gear may have a helix angle of zero to ninety degrees, theoretically. 
However, the range of these halix angles may be cut down by certain practical design con- 
siderations, such as pitch or pitch diameter size. For instance, a speedometer drive gear 
must have a hole big enough to fit over the transmission output shaft, which limits the gear 
pitch diameter to a minimum value and in turn limits its helix angle. 


The second part of this problem is to find the load-carrying capacity of the gears selected. 
Since the load-carrying capacity is a function of the allowable maximum stress induced by 
the load, it will be the object of this paper to show how both the maximum surface compres- 
sive stress and the deformation at the point of contact may be found on a spiral gear set for 
any given load. This particular analysis is applied to involute helicoids only and cannot be 
applied to enveloping worm gear sets or to Gleason-type hypoid gears. 


It should be noted that the theoretical instantaneous surface contact which occurs between 
crossed-axes gears is that of a point. In parallel-axes gears this theoretical contact is a 
straight line running the full face width of the gears. Because of this point contact, the load- 
carrying capacity of spiral gears as compared to parallel-axes gears is very low. Thus a 
comparatively small load produced a relatively high compressive stress at the point of con- 
tact. This point of contact deforms under load to a small ellipse of contact, and the two con- 
tacting gear teeth actually come closer together due to this flattening. On most gear sets 
with higher crossed axes angles the total deformation of the teeth that occurs is due solely to 
this surface flattening. The tooth bending that takes place is usually negligible due to the 
low tooth load. Likewise the bending stress occurring at the root of the gear teeth is small 
enough to be ignored. Any failure that will occur will be at the points of contact only, such as 
failure of the surface and not bending failure. 


Quite a few years ago, gear design experts suggested using the Hertz equations of contacting 
surfaces for determining the stresses and strains for parallel-axes gears. The work done by 
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Almen, Buckingham, McFarland, and Straub and the present A.G. M.A. work are all based on 
the Hertz formulas. This method of attack will be continued here on crossed-axes gears.! 


The validity of the Hertz equations requires several assumptions, which must also hold true 
when the equations are applied to spiral gears. These assumptions are: 


(1) The materials in contact must be homogeneous. 

(2) Contact must take place under static conditions. 

(3) The subsurface stresses in both bodies must decrease to zero before reaching the out- 
side diameter of a gear or the opposite tooth flank surface or either face. 

(4) The proportional or elastic limit of the materials must not be exceeded. 


The first condition of homogeneity and the fourth condition of limiting stress are easily met, 
because most metals are sufficiently uniform to be considered homogeneous (excepting per- 
haps some cast materials) and because we do not desire to go beyond the elastic limit as that 
is considered to be failure. For the third condition, it is only when the contact occurs near 
the outside diameter of a gear that the subsurface stresses do not become zero before reach- 
ing the outside diameter; otherwise most gear face widths and tooth thicknesses are suffi- 
ciently large to fulfill the requirement. However, for the second assumption it is impossible 
for running gears to be in static contact. Until someone develops other than an empirical 
answer we will have to use velocity correction factors, which at the present time are almost 
universally used on parallel-axes gears. Such are the PV or PVT maximum limits or the 
Barth type factors. 


Maximum and Minimum Radii of Curvature of Involute Helicoids 








It is a known fact that all involute helicoids can be generated by a series of straight lines that 
are all tangent to a helix on a base cylinder, as shown in Figs. 1 and 2. Hence it can be said 
that an involute helicoid is a surface composed of a series of straight lines called gener- 
atrices. For any given gear helix angle and normal pressure angle, these generatrices have 


Figure 1. “Involute Helicoid”. Straight 
Generatrices in Helicoidal Surface and 
Involute Profile of its Transverse Sec- 
tion. Courtesy of National Broach & 
Machine Co, 





1A first step in this direction was an approximation published by Prof. E, Buckingham on 


pe. 534 of Ref. 1. For details see our footnote 2. 
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a given angle to the gear axis. This angle is the base helix angle, which is the angle of the 
helix on the base diameter. The multiple relationships between the helix angle ¥Y, the normal 
pressure angle @,, and the base helix angle Yh can be found in Dr. W. F. Vogel’s “Involu- 
tometry and Trigonometry,” page 266 of Reference 8. 








Figure 2, Tangency of Straight Generatrices to 
Base Helix of Surface. Courtesy of National 
Broach & Machine Co, 


To use the Hertz equations, it is necessary to find both the minimum and maximum radius of 
curvature of each of the surfaces in contact and also the angle between the planes which con- 
tain the principal curvatures of both surfaces. 


We have already seen that the maximum radii of curvature of both members are infinite. 
Their arcs are contacting straight-line generatrices of the two involute helicoids. 


When two crossed-axes involute gears are in mesh their contacting generatrices, the straight 
lines, will intersect at some angle H, since the generatrices coincide only for parallel-axes 
gears. Therefore it will be necessary to find the angle H. 


The minimum radius of curvature is required next. We know from the theory of continuous 
surfaces that the minimum radius of curvature of a surface at a given or selected point Q 


lies in a plane which is at right angles to the plane containing the maximum radius of curva- 
ture of the surface at that point. 


Therefore we have a series of straight lines which compose the surface of the helicoid and 
another series of lines at right angles to this surface which will be the normals of the sur- 
face. These latter lines will be tangent to the base cylinder, and the angle they produce with 
the axis of the cylinder is the base lead angle, or 90° minus the base helix angle. The length 
of such a normal measured from the point of tangency on the base cylinder to the point Q on 
the surface is the minimum radius of curvature at the given or selected point Q. 


From Figure 3, it can be seen that the minimum radius of curvature of the surface is Pq at 
any point Q of the involute helicoid whose transverse pressure angle is Dq: 


_ R, tan ~. 
g Cow W 
¥, 
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NVOLUTE GENERATRIX 


HELICOID \ 


Pd Figure 3. Principal Radii of Curvature of 
’ Involute Helicoid. The maximum radius 
INVOLUTE YY / of curvature ofan involute helicoid (name- 


ly, that of any ofits straight generatrices) 
is infinite in magnitude. The minimum 
radius of curvature of this surface at the 
freely selected point Q is labelled Pa: By 
contrast, BQ is the radius of curvature 
of the involute at Q but not a “principal” 
radius of curvature of the helicoid, Note 
that the normal pressure attacks in the 
direction of Pq, if Q is the contact point of 


a helical involute tooth. 


Y” “BASE HELIX 
TRANSVERSE PLANE 


Radii of Curvature of Varying Contact Elements of Meshing Spiral Gear Tooth Surfaces 











For meshing gear teeth it is necessary to find Qq values of such surface elements which are 
instantaneously in contact. 


The values gq, and pq, of two contacting spiral gears may be considered as the radii of two 
steel cylinders, crossing each other at an angle H and pressed against each other with a force 
F,, equal to the “normal” tooth load. If angle H is established according to the equations of 
the next section, then the contact of the “substitute steel cylinders” is equivalent to the point 
contact of the spiral teeth for all calculations of stress and strain. ? 


It can be seen from any involute spur gear set that the total theoretical length of the path of 
contact is constant for any given operating pressure angle. Therefore, the sum of the radii of 
curvature of the gears at any point of contact is equal to this theoretical length of the path of 
contact. This will also hold true for crossed-axes gears. The sum of the two values of 
contacting surface elements will remain constant for any spiral gear set. Then as the Pq of 
one gear decreases, the py of the mating gear must increase. Therefore, it is necessary only 
to find the sum, and then it will be possible to find the stress or strain at any point on either 
gear profile. The equation for this sum is obtained by finding the two special pg values at the 
pitch point, which we will label p, and p,. If, by the same token, we label the known trans- 
verse main pressure angles at the pitch point @, and @,, then we obtain from Eq (1) 


2This concept of “Substitute Rolls” has been introduced onpg. 534 of Ref. 1. But the approximate 
haracter of this solutionleadsto errors of both the “Roll Size” andthe “Angle H” between the 
rolls. The error of the resulting tooth load is always on the safe side. 
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This (p, + p,) can be considered as the length of a piece of string which is tangent to both 
base cylinders when they are placed at their proper center distance and crossed axis angle. 
The string must also be tangent to both normal base helices. 1 * p,) is therefore equal to 
the distance from the tangency point on one base cylinder to the tangency point on the other 
base cylinder. 


The greatest Pq values occur at the outer diameter of each gear and will be designated by 
Pout- By the same token, the smallest pq values occur at the innermost contact radius and 
will be labelled pj,. They can be easily determined by the difference between poyt and 
(p, + p,). These special values pin and Poyt are limiting values of Pq and mark the beginning 
and end of gear contact. They may be found as follows: 


Vii - Bay 


= (3) 
Joust Cry, 





Jnl * (P, tf) Pose t (4) 


[5% 2 
p _ Roz - Be ts) 
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Prinz * (ft 2) fut 2 (6) 


In these equations Rp, is the outside radius of the gear and Rp is the radius of the base cylin- 
der. 


After finding the limiting poyt and pjn values, the stresses and deflections at the limiting 
contact positions can be computed. Then we can arbitrarily select as many additional values 
of Pq between these limiting values as we desire to compute the stress or deformation. 


Angle H Between Contacting Generatrices of Meshing Spiral Gear Teeth 





This angle H is constant, i.e. independent of the contact position of the gears in any gear set. 
Its equations, derived for three types of gear sets, are: 


Angle of gear crossed axes (shaft angle) 

n = Normal pressure angle ) 
Transverse pressure angle | at the pitch cylinder 
Helix angle 


<€oom 
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WY, = Base helix angle 


All other angles are parameters. 


Then for = 90° ¥, + ¥, (Right-angle drive) 
tan B tan», sin Q, 
tanH = tan (8 + @) sin Vp 

For x ¥, -¥,, whereY, >, (Spiral set having opposite hands) 
tan a tan Wp, cos Q, 
tan B tan Wp, sin @, cosa /sin (a - &) 
tan H tan (8 - @,) sin Vv be 

For x % (a helical and spur gear meshing together), where Y, = 0. 
tan H sin WY), tan O, 


For the most common occurrence, 90° drive sets, the chart (Fig.7) shows the angle H for any 
helix angle for 14-1/2°, 20° and 25° normal pressure angle gears. The use of this chart saves 
the computation of the above equations. 


The important aspect of this chart is that it is possible to get a yeneratrix intersecting angle 
as low as 28°, which is equivalent to having two “substitute cylinders” crossed at 28°. It is 
obvious that their load-carrying capacity is much greater than that of the “substitute cylin- 
ders” placed at 90°, as advocated in the approximate solution of Reference 1. 


Details of Derivation of Angle H 





The derivation of the equations for a right-angle drive gear set is shown in Fig.4. A is the 
true length of the generatrix of gear 1, whose projection in the transverse plane is inclined at 





Figure 4. Illustrating the Derivation fora Right~ 
Angle Drive. The involute helicoids of two 


meshing teeth contact each other in only one 
point of their helicoidal surfaces. The illustra- 
tion shows the establishment of the angle H be- 
tween the contacting generatrices of the two 
helicoids. 
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the transverse operating pressure angle @,. The angle p, is the base helix angle of this 


gear. The angles @, and ¥ b, are the corresponding angles onthe mate. Then from Fig. 4 


of the preceding page 
; : i 
ba V; hin 2 7 
Me F bi - 1 
Lan f= a OE 5 tam Hf Anh, (4) 
A cond, 





Then H = Ci 
Auat;, or(b+ $¢,) 
Cor /2 Ain Uy 





elias my ® 
Tan p A $, ) AAn 32 
or eliminating the parameter B, 


: (Can Vy pm p, + lam ?, jun Vor (8a) 
- Tan Y,, Aun P, Tan Pr 





Tam H 


The derivation of the angle H for the gear set with crossed axis angle equal to the difference 


of the gear helix angles is shown in Fig. 5. 


Tam Aan ty, md = an Fiy am $, ter a 
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ANGLE 2; 


Figure 5. Gear Axes Crossing at Any Angle. 
The establishment of angle H is shown for the 
case when the angle = between the gear axes is 
the difference of the helix angles 


and Y, at 
the pitch cylinders of the gears. 


1 





Jan He an La-Z)pin (4- b,)or Vs, 
ce 2 
Tan H = Tan (B-¢: ain Yr (11) 


The derivation of the angle H for a spur gear mating with a helical gear is shown in Fig.6. 





tam «= A din Vir OF = tan fy, ord, 
A ar Yj 
Tan #= A bin Vj ting, /02br _ Tan Vy Mng, 
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AR GEAR 2 
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% \/GEAR 1 (HELICAL) 
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Figure 6, Contact between Helical and Spur 
reeth. Establishment of the intersecting angle 
H of the contacting generatrices of the meshing 
tooth surfa< es. 


Zam Yj) 1 @; 
C0 o,, (APE + pr Dl lan Yj, @#,) 


Tan H = iin VY Zan D, (12) 





Maximum Compressive Stress and Strain 


The formulas, constants, and coefficients are taken from Hertz and later investigators. See 
References 2, 4, 6, 7 and 9. 


It should be noted that K is a constant influenced only by the materials of the gear teeth, while 
the other constants, 5 and @ depend on the selection of the contact point. See Eqs (13) to (15). 


© is an auxiliary angle, whose establishment yields the coefficients Cg, Ch and Cg by inter- 
polation from Table I. Additional explanations are given along with the consecutive steps (13) 
to (19) of the calculation: 





— A E, Ey 
er 7 es) (13) 
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Figure7. Chart for Right-Angle Involute Spiral 

Gear Sets. The curves give the intersecting 


angle H of the contacting generatrices of the 
meshing helical tooth surfaces. This angledoes 
























































not change with changing contact positions of 
the gears, 


where E is the modulus of elasticity of the material of either gear and p is the Poisson’s 
ratio of the material of either gear. 








Aane= 57 + 2anen (15) 
Oh "Apa 


where p, and p, are the various values Pout and pin taken from either Eqs (3) and (4) or Eqs 
(5) and (6) or any two Pq values belonging to any contact point between these limits. 


-C, ve 


}= 1 (AE 








TABLE I 


Coefficients for Equations (16) through (18). 
Excerpt from Table XIV of Reference 4. 


° Ca | Ch Cq 
0° —_ 0 aes 
10° 6.612 319 851 
20° 3.778 .408 1.220 
30° 2.731 493 1.453 
35° 2.397 530 1.550 
40° 2.136 567 1.637 
45° 1.926 604 1.709 
50° 1.754 641 1.772 
55° 1.611 | .678 1.828 
60° 1.486 117 1.875 
65° 1.378 .159 1.912 
70° 1.284 802 1.944 
15° 1.202 846 1.967 
80° 1.128 .893 1.985 
85° 1.061 944 1.996 
90° 1.000 | 1.000 2.000 


Where a and b are the respective major and minor semi-axes of the ellipse of contact, and P 
is the normal tooth load. 


-— 3 ee ) 
d=C, et (18 


where d equals the total deflection at the point of contact due to surface deformation. 


¢ 4s P 

re 19) 

C ( 
Tab 

where S,. is the maximum compressive stress. 


Thus by using the proper values of H and p it can be seen that it is possible to obtain the the- 
oretical stress and deformation in a crossed-axes gear set for any contact position of the 
gears. 


The problem of what is the maximum allowable surface stress can only be answered by tests, 
but the above equations allow the designer or tester to reduce the loads to a basic stress 
value. When these basic stress values are considered with the sliding velocity, it should be 
possible to get a formula for the maximum allowable load. 
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A GENERAL STEADY-STATE SOLUTION OF THE TWO-MASS TORSION PROBLEM 


FARL F, RIOPELLE 
Houdaille-Hershey Corporation 


Detroit, Michigan 


Introduction 


The two-mass torsional vibration system is widely used as an example of the application of 
mathematics to analyze the behavior of machines. Generally, existing articles and texts have 
analyses which give the natural frequency of the system, and the total relative motion of the 
two masses when only one mass is driven by an exciting force. 


Other factors can be obtained which are useful to the engineer who is concerned about the 
performance of such a machine that can be idealized as a two-mass vibratory system. These 
are: 


1. Relative amplitude of the masses when they are both excited by harmonic torques at an 
arbitrary phase angle. 

2. The absolute amplitude of the masses. 

3. (a) The nodal points of the system at a given forced frequency. 
(b) The nodal point at the resonant frequency. 


The mathematics required to derive expressions for the above conditions is elementary, but 
the results are useful. The derivation, results, and discussion are given in the article. 


Discussion 


This solution defines the torsional response of idealized driving and driven machines con- 
nected by a flexible means. Driving machines are exemplified by electric motors, hydrau- 
lic motors, gasoline and diesel engines, etc. Driven machines may be electric generators, 
fans, hydraulic pumps, air and refrigerator compressors, etc. The flexible connecting means 
are usually shafting, couplings, chain-gear-and-belt drives, etc. 


An analysis such as this is made to determine the loads or stress developed in the coupling 
medium by the torque fluctuations of the driving or driven machine. The vibratory loads are 
often neglected in the determination of the size requirements of the connecting medium with 
results of a seemingly unexplainable failure. 


The analysis is made by idealizing the system into two masses of polar moment of inertia 
equal to those of the driving and driven machines connected by a torsional spring of the same 
rate as the coupling medium. It is also necessary to know the harmonic content of the torque 
fluctuation of both the driving and driven members. 


The idealized system can be represented schematically by Fig. 1, in which J, and J, are polar 
inertias (inch pound second’), 6, and 6, are the angular displacements of J, and J, (radians), 
K is the torsional spring constant (inch pounds per radian), T, and T, are maximum torques 
of the harmonics considered (inch pounds), N and n are harmonic numbers, w is angular velo- 
city (radians per second),t is time (seconds), and a and # are phase angle relationships of the 
torques (radians). 
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The accuracy of the analysis is wholly dependent upon how accurately the above dimensions 
represent the actual machine. 
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Figure 1 
The next step in the analysis is to write the differential equations of motion of the idealized 


system. This is done by equating the sum of the external forces on each mass to the rate of 
change of momentum of the mass. Thus, 


46+ K(Q-O@2)= TF sinfnwtr«!) 


oF Se. +K(@,- 98)= Te sin (Nwt +3) 


“ _ Ji f 
OrF (8,~- 82) = F— sin (nwt +x) 
S.- *~/e -©@,) = fe. sin(Nwt +f) 
£ Je ‘ ™ Jez 
KX 4 
t 7 = 70° = s ; 
“a = 2 = fs 
Ty 7e 
7 alan Je = Qe 


Go-& + (ptt (2,- @2)= 9, sin nwt + x)- Q, sin (Nwt +) | 


The prime concern is of the relative motion of the two masses which can be used to determine 
the stress or load in the coupling medium. Therefore let 


a“ 


9-e,-$ aed 6-2,=¢ 
Thus. $+ (p+ 7h) P= eo) sin(/nwt+«<]}- OQ» sin (Nwt +) 
Let: fe + Hz = pP* 
g wx P* = QYsn (nwt +a) - GQ, sin(Nwt +3) 
Expand the sine functions of the sum of two angles with an identity. 


pr PY = QO, [sin nwt cosx + cos mwl sina | 


-Q2 [sin Nwt cos G+ cos Nwt sin 3 |. 
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a and 6 are constants; therefore let 


cosec=-=A SiInN~= B 
cosfgf=C sin B=D 


p+ P*o QA S/*7 nwt a Bocos nwt]-Q[C sin Nwt +DOcos Nwt | 


The steady-state condition of the resulting motion is described by the particular integral of 
the above differential equation. The particular integral is easily found by the method given in 
“Mathematics of Modern Engineering,” Doherty and Keller, Vol. 1, pg. 41., equations 90, and 
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This is a convenient form for the expression of @. 


Notice when P= nw or Nw, the angle @ becomes infinite, which would cause failure of the 
coupling means. 


The condition of P = nw is known as resonance or the critical frequency and, when expressed 
in the expanded form, is 
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are known as dynamic magnification factors, and they can be plotted as a general curve 
Fig. 2). 


Obviously because this analysis describes the twist resulting from the torque fluctuations, the 
twist resulting from the uniform load must be added. 


0-0, +O 
© = total twist, @, = uniform twist, and @ = fluctuating twist. 


, nw Nw 
Observe that in order to reduce the effect of fluctuating torque, >> must be greater 


than V2. 


It is customary to design the coupling medium with a torsional (K) low enough to insure a dy- 
namic magnification factor of less than 1 for the principle harmonic of lowest number, thus 
insuring that all higher harmonics are automatically taken care of. 
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Thusfar no mention has been made of the effect of phase angles a and B. It is apparent that if 


sin (nwt + a) and sin (Nwt + B) equal 1 and -1 simultaneously @ will be greater than if they do 
not. 


It would be desirable to select a and B so as to minimize @ The surest way of making this 
selection is to plot separate sine waves of sin (nwt) and sin (Nwt) and plot their sum for a 
sufficient number of phase relationships between them to establish the minimum. 


Although the total twist ® is of prime concern there are conditions when the amplitudes, 8, and 
6, of the masses J, and J, are necessary. This would be particularly true if J, or J, werea 
reciprocating engine or compressor where the vibratory amplitude of the shaft would contri- 
bute greatly to the bearing loads on the unit. 


8, and 8, can be found from the information already obtained. The computations will be sim- 
plified if it is noted that in the expression for @ the effect of T,, T,, a, and B are independent; 
hence, they can be treated by the method of superposing their effects. Thus the resultant am- 
plitudes 6, and 6, as developed by T, will be shown; the same can be done for T2, and the re- 
sults added to the effect of T,. 
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Substituting, 
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Conditions for determining C, and C,: 
6, is a harmonic motion and goes to zero, thus, C, = C, = O. 


: 3m 5n 
8, is a maximum when sin nwt =1; therefore, nwt o> 


. and so forth. 
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Solve for 8, in a like manner. 
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Conditions for determining C, and C,: 
C,=0O 


6, is a harmonic motion, thus passes thru zero; therefore, C, - 
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8, is a maximum when sin nwt = 1; therefore, nwt >> and so forth. 
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It will be noted that the displacements 8, and 8, are in opposite directions when the exciting 
frequency is greater than the resonant frequency. This implies that there is a point in the 
connecting medium where there is no motion or displacement. This point is called the node. 
It is obvious that if a gear train or chain drive is to be installed on the connecting shaft it 
would be desirable to locate the drive at the node to eliminate back lash rattle which results 
in tooth fatigue. 


The position of the node can be determined by using the fact that the sum of the torques in the 
shaft is zero. As 


-K, ®-=K, &-K@ 


where K, is that portion of K from J, to the node(in lb./rad.),and K, is that portion of K from 
J, to the node(in lb./rad.). 


1 1 , 
Note that K = K, +£e (Fig. 3). 
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Reduce to lowest terms 
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It is an interesting revelation that if J, is excited by a frequency (nw), the node is formed at 


the position that would form a single-degree-freedom system for J, of natural frequency (nw). 
See Fig. 4. 


This makes it very simple to determine the position of the node for any exciting frequency. 


The node at the resonant frequency of the system will be as follows: 
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Note the resonant frequency is K a . 
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TRANSIENT MOTIONS IN LINEAR SYSTEMS 


JT, ORMONDROYD 
University of Michigan 
Ann Arbor, Michigan 


Most readers of this lecture are familiar with simple linear differential equations, such as 
mxX + cx + kx = f(t) 

either in this form or in the old operational form 
(mD°+ cD+ k)x = f(t) 

More complicated equations of the same nature exist which can be written 


Z(D) x = f(t) 
a-l a-2 
D DY 


There are three aspects to the solutions of equations of this type. 








where Z(D)=(... a,D’+a,D+ a+ 


(1) Z(D) and f(t) may be given completely and explicitly, the problem being to integrate the 
differential equation so that X(t) is found as an explicit function of time. This might be 
called the direct problem. 

(2) Z(D) and X(t) may be given, the problem being to infer f(t). This might be called the 
inverse problem. 

(3) X(t) and f(t) may be given, the problem being to find Z(D). This is the problem inherent 
in Carson’s Integral. 


Problem 1 is the classic problem of all textbooks and the LaPlace Transform Method is un- 
doubtedly the most acceptable way to carry out the solution for X(t). 


Problem 2 is seldom discussed in detail in the textbooks. One common physical problem 
which leads to problem 2 is the basic problem of all instrumentation in which the mechanical 
and electrical circuits are all linear. The fundamental problem of all such recording in- 
struments is: 


Given the characteristics of the instrument, Z(D), and the record X(t), which might be desig- 
nated the output of the instrument, what is the correct and accurate value of the disturbing or 
input function, f(t), the quantity that is being measured? 


Problem 3 arises in a calibration of instruments mechanically and electrically so complex 
that no complete analysis of the characteristics of the instrument exists. This type of prob- 
lem is practically never discussed in detail, although Carson’s Integral is the expression of 
its formal solution. Carson’s Integral holds when f(t) is a step function and X(t) is the re- 
sponse A(t) to the step function. 
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Problem 1: The Direct Problem 





While the LaPlace Transform Method is the best method for the solution of this problem, 
it is frequently useful to use the older Duhamel Integral Method on some problems. This 
method is discussed by Timoshenko in Vibration Problems in Engineering, 2nd Edition, pp. 
102-109, and by Karman and Biot in Mathematical Methods in Engineering, Chapter X. The 





emphasis in these treatments is put on systems of one degree of freedom. But the Duhamel’s 
Integral method is even more useful in systems of many degrees of freedom. 


If an impulse of short duration is applied at any point in a mechanical system, that point ac- 
quires a velocity which is added to the velocity already existing at that point. The added ve- 
locity is independent of the configuration of the system at the instant of application of the im- 
pulse. The added velocity is 


At 


Vv Fdt 


Ble 


0 


at 
in which m is the mass of that section of the system in which the impulse J Fdt is applied. 


A unit impulse is defined as - Fdt = 1 lb-sec as F—» om and at—»o. The velocity caused 
by a unit impulse is V = 1/m. The response of the system to a unit impulse is designated as 
hj,t,j,7 in a lumped system and as hx t.c7 in a lineal distributed system. In words, hit jz is 
the displacement of body i of the system at time t caused by the application of a unit impulse 
to body j of the system at time t= 7. For the lineal distributed system, hx t ¢,7 is the dis- 


placement at point x and time t caused by the application of a unit impulse at point x = c and 
at time t =T. 


To illustrate the use of Duhamel’s Integrals, let us take two systems shown in Figure 1 and 
Figure 2. 


Lumped System 





The equations of motion for the system shown in Figure 1 are: 


LLLLLLLLLLLLAL 








=k, 
m,x, + (k, + k,)x, - kx, =O 
m, X, 
K -k, x, + m,X, + (k, + k,) x, - k,x,= f(t) 
m, * 
-k, x m,X, + k,x, =O 
f(t) SA, , 
Mm, Xs 
Figure | 


To apply Duhamel’s Integrals, we solve the easier problem in which f(t) = O for the initial 
conditions at t = 7 of x,= x,= x,= O, X, = x, = O, and x, = 1/m,. 


These initial conditions will arise if the system has a unit impulse applied to body 2 at times 
t=y7. The system is otherwise assumed completely quiescent at that time. 
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The differential impulse applied to body 2 is f(t)dt. The actual responses of the bodies in 
the system to the exciting force f(t) are: 
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These quadratures give the results of a system initially at rest when the disturbing force f(t) 
is applied. This case occurs frequently in physical cases of transient motions. 


Distributed Lineal System 





Considering bending only, the equation of free vibration for a simply supported beam of uni- 
form cross section is 
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Figure 2 














x,t = displacement of section at x at time t 
Young’s Modulus of beam material 
weight per unit volume of beam material 
rectangular moment of inertia of beam cross section 
area of beam cross section 
= acceleration of gravity. 
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If at time t = O, for all values of x, y= O, then 
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where 


If at time t = O a finite velocity V is applied to small length 6 on the beam centered at x = c, 
this condition can be shown in Figure 3. 
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By the use of the ordinary Fourier Integral, 
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If the initial velocity at x = c was caused by a unit impulse, 


as 
ea 





and 





where m = entire mass of the beam. 


For a unit impulse applied at x =c and t =f, 
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If the problem concerns the deflection of the side-rod of a locomotive, the inertia load on the 
side-rod is a function of position and time W = Wosinwt _ Ilb/in. 


where Wo is the maximum value of the distributed load. 


The differential impulse is W,dc sin w} dy. Note that sinnnc/f is the “influence function” 
for this problem and the variable of integration is c, not x. 
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Problem 2: The Inverse Problem 





Problem 2 is discussed at some length in the article on “Motion Measurement,” Sections 40- 
43, pp. 375-380, in the Handbook of Experimental Stress Analysis. It can be illustrated by the 
simple seismic instrument shown in Figure 4. 
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The equation of motion for this system is 


n(4+X)+ex+kxe= O 
or 
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Figure 4 
where y(t) the motion being measured 
x(t) the relative displacement between mass and case 


damping constant 
critical damping constant 





spring constant 
mass 


BIA O|O 


Seismic instruments respond only to accelerations, and the output is only the relative motion 
between mass and case. If the relative motion x is magnified mechanically, electrically, or 
optically and transferred to the record in a linear fashion, the record x is 


X = px 


Therefore the record X is related to the input y as follows: 


YtenpX+pXs eg 


“ae Je 
A. 2X OX SF 


or . 





The record X os accelerometers is usually considered by the users of the instrument to be 
proportional to the measured acceleration. But the equation above indicates that this is only 
true if X and X are smallor P is very large. If none of these conditions exist, the record 
must be corrected by the addition of the terms x/p* and 2n/p X. These terms call for differ- 
entiating the record once and twice, instructions easier to give than to carry through. 


The above equation illustrates a common dilemma in all motion-measuring instrumentation. 
As p ~o the record X becomes more and more accurately proportional to y, but the sensi- 
tivity approaches zero. That is, the record becomes very small. Accuracy and sensitivity 
are mutually antagonistic goals. 

Suppose that the transient motion being measured starts at the beginning of the record so that 


a pega Xe X20 


Then the record equation can be integrated term by term and 


X + rep [xore er f[[xmdt= Hy 
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The record X is proportional to y only if f Xdt, [f Xdt dt, and p are very small. If none of 
these conditions exist, the record X must be corrected by the addition of the two terms 
2npjXdt and p’ //Xdtdt. These terms call for the integration of the record once and twice, 
instructions which can be carried out. All displacement meters have small values of p. 


Velocity Meter with Low Impedance Output 





A seismic instrument widely used is the velocity meter with an induction pickup coil and with 
recording made by a low resistance electro-magnetic oscillograph. This instrument consists 
as a whole of two mechanical systems connected by an electrical system. Each mechanical 
device can be approximated by a single degree of freedom system, and the electrical circuit 
can be approximated by a single or double mesh circuit. These are shown in Figure 5. 
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Figure 5 


The equations for this system are: 


n 


1) MX+CX+Kx+ fA, =- my 
2) 1, , +R, &,+ 2{(a,-as) At = dx 


3) 4,a, +R, a, - EA, aht ape o =O 


b) T p+Agr+deg=-VAz 


In which y = motion being measured, x = relative motion between mass and case in velocity 
pickup, @ = angle of rotation of oscillograph mirror, a, = current in mesh that is in the ve- 
locity pickup, a, = current in mesh that is in the oscillograph, m= mass of suspended weight 
in the velocity pick up (also a magnet), C = mechanical linear damping constant in velocity 
pickup, K = linear spring constant of weight suspension in velocity pickup, I = moment of in- 
ertia of mirror, in oscillograph, B = mechanical torsional damping constant in oscillograph, 
5 = torsional spring constant of mirror suspension in oscillograph, L, = inductance in velocity 
pickup and half the transmission line, R, = resistance in velocity pickup and half the trans- 
mission line, L, = inductance in oscillograph and half the transmission line, R, = resistance 
in oscillograph and half the transmission line, C = total capacity in velocity pickup, oscillo- 
graph and transmission line, € force per unit current reacting on suspended mass of ve- 
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locity pickup, ), = volts per unit relative velocity generated in velocity pickup coil, y = back 
electromotive force per unit angular velocity in oscillograph, and v = torque per unit current 
in mirror suspension of oscillograph. 


If p is the optical lever length in the oscillograph, the record is p @. The differential equation 
relating the record, the constants of the entire system, and the velocity being measured is 
(provided @ and y and all their time derivatives were zero at the time the record starts) 
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and the T, are characteristic times of the system. 


This equation indicates that the record may be a very much distorted picture of the velocity. 
The user of a velocity meter directly connected to an electromagnetic oscillograph merely 
hopes that every term on the left hand side is small except 1, nd In order to check the ac- 
curacy of the instrument, the user is instructed by the differential equation to integrate the 
record three times and differentiate it five times. He is instructed to know the characteris- 
tics of the instruments and the circuits involved in them. The integrations are possible. The 
differentiations are not. In practice none of these things are done. The operator hopes that 
the important components of his record will have frequencies lying between the low natural 
frequency of the velocity pickup and very much higher natural frequency of his oscillograph 
element. For the rest he trusts to luck. 


Problem 3: Experimental Method of Determining the Constants of an Unknown Instrument. 








In general, any seismic instrument which contains only linear mechanical and electrical ele- 
ments has the differential equation 


(rave arenretze BE eo)Xery 


or Z(D)X=-4 


If the analysis of the instrument and its circuits is too difficult to carry out, the values of the 
an can be found by an experiment proposed by Dr. R. D. Specht. This procedure depends on 
the evaluation of Carson’s Integral 


Qo 
Dz(d) if eA) dt 


where A(t) is the response of the instrument (the record) to a unit acceleration, applied to the 
instrument. If we are satisfied with unit accelerations less the g (the acceleration of gravity), 
it is easy to imagine and to approximate in construction a setup which will give constant 
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accelerations. For instance, the instrument might be mounted on a carriage which is accel- 
erated by a falling weight. If the instrument (pickup, amplifier, transmission system, and 
recorder) has sufficient damping, the response A(t) will settle down to a constant value, 9, 
within a finite length of time, T. There are only two possibilities, either @ = O, or 6 # O. 
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It is interesting to note that when © = O no integral terms exist beyond the first. Referring 


back to the discussion on the velocity meter, the equation relating the record to the disturb- 
ing function for accelerations is 


4 ~ 3 
(D+ - fe oie + , | ‘oan bb — an a 
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Now, a velocity meter subjected to a step acceleration is so constructed that @ = O. Here we 
seem to have a contradiction between Z(D) calculated by detailed analysis and Z(D) got from 
Carson’s Integral operating on a record taken with a step acceleration excitation. But if the 
detailed analysis is carried out on the velocity meter, it will be found that 
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so that there is no contradiction between the two developments. 
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SOME COMPUTING PROBLEMS IN THE AUTOMOTIVE INDUSTRY 


PAUL T. NIMS 
Chrysler Corporation 
Detroit, Michigan 


This first meeting of the Association for Computing Machinery in the Detroit area affords 
a welcome opportunity for some of us in local manufacturing industries to learn what is 
being done in the computing field and what kinds of machines may become available in the 
future. It also gives us a forum in which to discuss industrial mathematics and its com- 
puting problems. 


I have been asked this morning to introduce some typical problems of the automotive segment 
of Detroit industry. The term “automotive industry” is rather a broad one and means many 
things to many people; not only the vehicles that you and I drive on Sunday afternoons but 
trucks of all sizes, tractors, tanks and self-propelled guns, railroad equipment, internal- 
combustion engines for aircraft service, and perhaps even certain heating, air conditioning, 
and refrigerating equipment. No one man can be acquainted with all aspects of such a diverse 
collection of automotive products, to say nothing of the other things made in Detroit. Hence, 
this morning I shall mention only those matters that I have heard about and those few with 
which I have had some direct experience. 


This is an opportune time to talk about computation in the automotive industry for a number 
of reasons. There is the matter of the balance between experimental investigation and com- 
putational investigation. Where the physics of the problem is well enough understood to 
permit calculation of the strength and performance of the device with reasonable confidence, 
then the decision as to whether to calculate these matters for a range of designs or whether 
to test a series of designs is largely an economic one. As the cost of experimental proce- 
dures and of experimental parts is growing daily, there is a growing realization of the ad- 
vantages to be gained in time and money by substituting the computational approach for a 
large part of the usual experimental testing. 


There is also a special reason for talking about this kind of thing in the automotive industry 
instead of, let us say, the building industry. In an activity where it is traditional to make 
extensive changes in the product annually and where in the lifetime of one individual several 
sets of drastic changes are likely to occur, the burden of engineering work in connection with 
these changes becomes very heavy. For example, the introduction of a completely new auto- 
motive engine together with a new transmission, suspension, and braking system represents a 
whole collection of new engineering problems that must be solved rapidly and repeatedly until 
an optimum design of each of the components and of their interaction is found. 


There is still another reason why computation is important in the automotive industry today, 
and that is the flood of military contracts which we are beginning to face. While most of 
these devices will have rotating and reciprocating mechanical parts, the similarity to the 
peace time design problem ends at this point as these military devices are subjected to quite 
different stresses, weight limitations, and performance requirements. They also are often 
much larger. This means the design of completely new and often quite unfamiliar devices. 


For these reasons, therefore, I think the computing problems of the automotive industry 
are of special interest at this time and will mention a few of them both by physical kind and 
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mathematical kind in the hope that you will find the time to look into them further and perhaps 
devise equipment and methods which can be of general use in the industry. 


Kinds of Problems 


A. Substitution in Formulae 


Probably the most common recurrent problem is that of putting numbers into formulae. 
These formulae may relate to the performance of the vehicle, to the stresses in various parts 
of the structure, or to the necessary size of some parts. A man designing the differential 
gear train in the rear axle may cover pages with numbers relating to the shape of the teeth 
and the setting of the various knobs and dials on the machinery necessary to cut these teeth. 
The formulae abound in sines, cosines, arc tangents, and various products and quotients. The 
calculations for a single formate hypoid gear require the evaluation of about 670 items of 
which roughly one-third are trigonometric functions. The data are calculated to six figures. 
This calculation, with a desk machine, requires the full attention of a man for six days and is 
usually checked by a duplicate calculation by another man. While the programming of this 
problem will be difficult, the time saved should be considerable. 


The gear designer is also interested in smaller jobs of calculation. After agear is designed, 
it is often necessary to check the design for involute interference and for involute under- 
cutting. In addition, some special dimensions are required for the inspection of the gear 
after it is cut, a typical calculation being the diameter over inspection rolls. You may find a 
portion of the formulae which are used in gear design in a book entitled “Involutometry” by 
Professor Vogel at Wayne University. 


The design of the cams which open the valves in a modern automotive engine is an exercise 
in the juggling of acceleration forces and has been attacked with many quite complicated 
formulae. The designer is given the required lift, duration of opening, and permissible ac- 
celeration and must produce a table showing the lift to the nearest one hundred-thousandth of 
an inch for each degree of camshaft rotation. This table is:then copied on the drawing as a 
guide to the machine shop in the production of a master cam, and the copy must be rechecked 
for errors. A cam design method mentioned by Professor Winston Dudley in the January, 1948, 
SAE Transactions uses polynomials to represent the cam contour. This type of formula has 
been evaluated on punched card machines by several groups of people and in at least one case 
involved summing terms up to the twentieth power. 


Another very common calculation is that for finding the loads on the bearings of the crank- 
shaft. In this case, one must consider the magnitude of the force of the expanding hot gases 
as the engine goes through its thermodynamic cycle, the inertia of the piston and other re- 
ciprocating parts, the inertia of the rotating balance weights on the crankshaft, the angularity 
of the connecting rod (which gives a nonsinusoidal motion to the piston), and may even have to 
consider the deflections of the crankshaft. Such calculations must be made for an assortment 
of crankshaft angles covering two full engine revolutions, and all the forces must of course be 
summed with proper regard for their vector properties. This sort of thing can be done, and 
is being done, with an adequate supply of blank paper and a slide rule, but in view of the fact 
that it would be desirable to study these bearing loads at a variety of engine speeds and throttle 
openings, the computing problem assumes large proportions and these tools become inadequate. 


Some portions of the design of hydrodynamic torque converters are extended exercises in 
analytic geometry. A part of the specifications on the drawing is a table of the coordinates 
of the points of intersection of the ruled surfaces of the blades with the parabolic cone surface 
of the hub. The repetitive arithmetic of this calculation should be wel! suited to an automati- 
cally sequenced calculator. 
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B. Sequential and Iterative Calculations 


For some kinds of calculations we are not so fortunate as to have explicit formulae for the 
quantities desired and find it necessary to calculate through a set of formulae several times 
in order to find a satisfactory set of values for the unknown quantities. The design of an 
ordinary round wire spring is a simple example of this sort of thing, as it is necessary to 
make compromises between allowable stress, space occupied by the spring, the rate of the 
spring, its solid height, and the size of the wire used in its manufacture. Presumably the cri- 
teria for the compromise could be inserted in a setof instructions and the formulae evaluated 
repeatedly until a satisfactory solution is found, i. e., until all of the quantities desired are 
smaller than pre-assigned upper limits or until the machine stops and reports that no such 
spring can be built in the space available. 


The gear designer again seems to get the palm for complication when he attempts to modify 
gear tooth forms to obtain an optimum balance of tooth life between the gear and the pinion 
ina train. Present practice is to deal with the transcendental equations involved by various 
graphical methods which are quite time consuming. Programming of this calculation for an 
iterative solution method will be difficult, but easy computation may permit a more economical 
use of the alloy steels required for gears. 


C. Table Look-Up 


As a part of substituting numbers in formulae, one encounters over and over again tabulated 
functions of one sort or another. Trigonometric functions are quite common in any problem 
involving geometrical concepts, such as those just mentioned, and logarithmic functions 
appear frequently in the study of thermodynamic cycles. In brake drum analysis, thermal 
transients can best be described in terms of the error function so that tables of this and of the 
exponential are sometimes of interest. Now, of course, looking up numbers in tables is not 
exactly heavy manual labor, but it is possible for the volume of the work to build upso that the 
dangers of errors of transcription are great. These operations of consulting a table occur 
generally in the midst of the evaluation of formulae so that they are unusually awkward for 
machine calculation. 


D. Data Reduction and Automatic Data Taking 


Another class of substitution in formulae arises from the experimental side of the industry. 
When testing engines, for example, it is customary to record the torque, speed, number of 
seconds to burn a pound of fuel, and the temperature and barometric conditions. From these 
data can be calculated the corrected horsepower and the corrected brake specific fuel con- 
sumption together with other similar data. The amount of this kind of work varies with the 
kind of engine testing program which is under way and seems to be at a maximum with large 
aircraft engines. It is possible that some consideration should be given to the automatic re- 
cording of these data directly from the motions of the various indicating instruments and in 
such a form that the calculations can also be carried through automatically. This problem is 
complicated by the fact that atypical installation may have fifty different test cells in operation, 
and it would be desirable to record and calculate these data for any or all of the test cells 
without the use of separate, duplicate equipment in each room. 


E. Sorting 


Some of the problems of the automotive industry are quite similar to those encountered by 
chemists. For example, keeping a cross index of the system of internal reports by which 
the work of the engineering laboratories is recorded and keeping bibliographies generally is 
related to the problem of hunting through infrared spectrograms in an attempt to find spectral 
peaks corresponding to those of an unknown compound. This is primarily a sorting problem, 
although the application of the principles of logical truth calculators might make for much 
more intelligent selection of the bibliographical items to be given detailed study. 
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F. Sets of Algebraic Equations 


Automotive structures are generally indeterminate. That is, a V8 crankshaft is supported 
on five main bearings, and a modern steel body structure, while it bears some resemblance to 
a bridge truss, is very definitely not pin connected but typically has additional bracing at all 
corners. Study of the stresses and deflections of structures of this kind leads to sets of si- 
multaneous algebraic equations. Similarly, spectroscopic analysis of the components of an 
antifreeze solution yields a simultaneous set of algebraic equations, and rapid methods of 
solving these sets are convenient in speeding up the production of useful answers. 


G. Integration 


Thus far, the mathematical procedures | have mentioned are not beyond the attainments of 
a college freshman. This was deliberate. While occasionally it is necessary to make in- 
vestigations of the greatest mathematical sophistication, the burden of computation is largely 
that of lifting numbers and putting them into cubby holes, that is, of substituting in formulae. 
However, such matters as centrifugal stresses in rotating discs may require step-by-step 
integration for discs of arbitrary cross section, and test records of acceleration versus time 
may require similar treatment to find the velocities. The ordinary differential equation can 
appear in a study of the motion of an automobile and its passengers on a rough road. In fact, 
very complex sets of nonlinear differential equations will be required to discuss adequately 
the effects of shock absorbers, friction between the leaves of springs, and the special kind of 
spring rates which are sometimes used in suspensions. 


Extensive studies of the torsional vibrations of crankshafts and similar structures are made 
by some companies, and one of the members of the Industrial Mathematics Society reports 
that most of his professional activity is concerned with the repeated solution of the problem 
of two flywheels connected by a flexible shaft. Similar sets of ordinary differential equations 
may be obtained in the study of the transients in automotive ignition equipment. 


The partial differential equation turns up from time to time, mostly in the study of heat flow 
problems, such as the thermal conditions around a spark plug, which lead to a three-dimen- 
sional diffusion equation. There are also some possibilities of this kind in the study of aero- 
dynamics of grilles and radiator ornaments. 


I think it may be seen from this partial list that the problems of interest in Detroit are some- 
what different from those which inspired the first flowering of machine computation in this 
country. That is, we have relatively few occasions to solve a simple set of differential 
equations for a very large number of initial conditions. However, the problems are more 
complicated than that of computing payrolls so that we see here an intermediate computing 
problem between those already so well provided for in existing machinery. 


Machines 


It would be profitless to tell this group how to design computing machinery, but it seems 
worthwhile to mention some of the desirable properties of machines for industrial computation. 
The differential equations of interest usuaily represent the oscillations of some mechanical 
system with many degrees of freedom, and the result desired is usually an estimate of the 
amplitudes and periods of the motion. This kind of problem is well suited to the capabilities 
of analogue machines like the Differential Analyzer at Wayne University and the various elec- 
tronic differential analyzers which are commercially available. 


The remaining problems seldom involve the repeated evaluation of small steps along a curve 
but rather the evaluation of a few more complicated formulae. This kind of work can be done 
very comfortably without the blinding speed of some of the newer computers. Much more 
important than speed is ease of application, that is, simple programming and rapid changes 
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from one problem set up to another. For the optimum use of machine computation by auto- 
motive engineers, the machines ought tobe directly accessible to the man who has the problem, 
and it should not be necessary for him to have long, specialized training in tricks of pro- 
gramming. This probably means that a large number of subroutines should be available by a 
simple code, not only to extract square roots, but to solve right triangles, interpolate in tables, 
sum series, and perhaps do often repeated groups of operations used in a particular problem. 


Rather than burden the computing machine itself with these many special program orders, 
it may be desirable to build a translating machine, that is, one into which a formula for eval- 
uation may be entered symbolically and in which the various parts of the formula are divided 
into arithmetic steps and the appropriate holes punched in a program tape for subsequent use 
in the computer. If we think of each mathematical symbol in a formula as a logical proposition, 
it is entirely possible that chains of relays can be wired to select appropriate groups of pro- 
gram steps from a tape or card storage unit. These groups could be combined by the coding 
machine in the order dictated by the formula. 


Convenient operation will probably also require floating decimal operation and storage ca- 
pacity (for a hundred 12-digit numbers and their signs) which can be readily changed to stor- 
age of two hundred 6-digit numbers with sign. An additional convenience would be some 
method of typing the name of a quantity on the answer sheet alongside the magnitude of the 
quantity, as often the answers do not come out as columns of pairs of values of x and y co- 
ordinates. 


Conclusion 


I have indicated some of the computing problems of the automotive industry and some de- 
sirable properties of machines for their solution. It should be emphasized again that the com- 
putations are probably too diverse to allow the use of special-purpose machines and that ease 
of application is quite important. Finally, if such equipment can be made available at less 
than the price of a new car, it should be of wide interest. 
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A SURVEY OF INDUSTRIAL SUBSIDY OF THE STUDY 
OF APPLIED MATHEMATICS 


THOMAS H, SOUTHARD 
Wayne University 


Detroit, Michigan 


In the spring of 1950, the Education Committee of the Industrial Mathematics Society wrote to 
twenty-four organizations representing various industries (automobiles, steel, oil, electrical 
manufacturing, utilities, chemicals, business machines, aircraft, railroads, radio manufac- 
turing, etc.) in an attempt to obtain information of use in the Society’s campaign to promote 
the use and point up the relative efficiency of the analytical approach to many of the problems 
of production, research and development and related industrial fields. Twenty replies were 
received. 


The twenty replies included four from the automotive industry, three each from oil companies 
and utilities, two each from electrical manufacturers, business machine companies, distillers 
and one each from the steel, chemical, aviation, and radio manufacturing industries. 


Our attempt to secure information included the following two specific questions: 


(1) What courses in “applied mathematics” (beyond the calculus, but including review 
calculus as well as such courses as hydrodynamics, elasticity, etc., not normally called math- 
ematics) are you offering for your employees on your own premises or elsewhere? 


(2) What subsidy, in the form of time off, money, or other forms of encouragement, do 
you give your employees while attending courses sponsored by your organization or offered 
by nearby universities? 


Seventeen of the twenty replies indicated some kind of subsidy of the study of applied mathe- 
matics. These seventeen affirmative answers indicate essentially two types of programs 
which can be characterized either as “regular” or “occasional”. 


Seven of the seventeen affirmative respondents are companies which sponsor “regular” (i. e. 
regularly scheduled) programs including work in applied mathematics. Four of the seven 
regular programs are degree programs, three are non-degree. 


The typical degree program leads, on completion, to a full-fledged college or university 
degree in engineering, mathematics, or physics. In one case an employee can work for the 
bachelor’s degree, in another the master’s, and two companies sponsor programs leading to 
the Ph. D. Most of these degree programs are staffed partly by university personnel, partly 
by company personnel. The typical program is in operation during the evening; some courses 
are taught on company premises, others at the cooperating university. 


Universities co-sponsoring these degree programs include Illinois Institute of Technology, 
University of Pittsburgh, Indiana University, New York University, Polytechnical Institute of 
Brooklyn, Stevens Institute of Technology, University of Pennsylvania, Cornell University, 
and John Hopkins University. Employees in some instances pay only part tuition, with pre- 
miums for good grades in some cases: in other instances they enjoy free or reduced tuition 
for courses taken under an Industrial Lecturer (company employee on the teaching staff) but 
pay full tuition for courses taken under the resident university staff. 
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In addition to the four companies with regular programs leading to a degree, three other 
companies have regularly scheduled non-degree programs. These are usually staffed by 
company personnel and are taught during the day on company time at no student expense. 
These programs are in general more likely to be aimed at training company personnel in 


techniques of special use to employees in their work, but they do include some applied math- 
ematics. 


Besides the regular programs, there are occasional courses sponsored by fourteen of the 
twenty respondents (this number includes four companies which also have regular programs). 
There are three typical kinds of subsidy involved: 


(1) A course on the premises, for employees, taught by employees or nearby university 
personnel, during working hours. 


(2) A course on the premises, for employees, taught by employees or nearby university 
personnel, after working hours. 


(3) Regular university courses, chiefly evening. 


For the first category the only subsidy is time-off (and the instructor’s salary when an 
outsider is involved). In the second case, the only subsidy is the instructor’s salary, if any. 
In the third case, the typical subsidy is part of the tuition, although in some instances even 
more is subsidized. This part tuition varies from 1/3 to 2/3 and in some cases the greater 
subsidy is made only for A or B grades. The courses must be approved by the company and 
must be rather directly related to work the employee is doing or will be doing, unless they 
are part of a degree program which satisfies this requirement. 


One company allows its employees up to five hours per week of time off to attend approved 
daytime courses at nearby universities. Another company occasionally sends one of its re- 
search men to summer school for more work in applied mathematics. 


As to specific courses mentioned by the several companies, the most popular appear to be 
Mathematics (or Calculus) Review, Ordinary Differential Equations, Advanced Engineering 
Mathematics, and Probability and Statistics. Not far behind are Advanced Calculus, Theory of 
Complex Variable, Vector Analysis, and Operational Mathematics. 


Most of the companies replyingto the questionnaire are enthusiastic about the results of their 
efforts in this area of their over-all training program. They feel that the subsidy pays for 
itself both directly and indirectly. 


Those who have gone furthest with their programs are inclined to be the most enthusiastic 


about the results to be attained. One letter ends as follows: “I ... believe that the program 
upon which your Society has embarked is well worth the time and effort expended.” 
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THERMAL DESIGN OF EVAPORATIVE CONDENSERS 


HUGH J. SCULLEN 
Consulting Engineer 
Detroit, Michigan 


The purpose of this paper is to show how fairly elaborate methods of analysis and synthesis 
may very profitably be applied to test data of a relatively low degree of accuracy when the 


object is a comprehensive rather than a precise picture of the character of a heat transfer 
device. 


The particular kind of device used as an example is what is known as an “Evaporative Con- 
denser”. This is essentially a forced-draft, air-cooled condenser provided with means for 
maintaining a layer of water on the air side of its heat transfer surface so that heat will be 
removed from it by evaporation as well as by convection. 


The specific apparatus used to obtain the test data used in this paper is shown in Figure 1. 
Details of the condenser are given in Table 1 and the actual test data are shown in Table 2. 
The tests were runas partof a preliminary study and, consequently, were run with the emphasis 
on speed with moderate accuracy rather than extreme precision. For this reason as well as 
because the data is used here simply to demonstrate a method of handling a heat transfer 
problem, no detailed description of the test set-up or procedure will be given. 
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SKETCH OF TEST SET-UP 





Figure | 


The discussion will be limited to a study of the design factors of heat transfer only. No 


mention will be made of such other very important factors as the air fan or water pump and 
their accessories. 
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The heat dissipated by an ordinary dry-surface, air-cooled condenser can be computed by the 
well-known equation 


q=u a 4u- dt fs (1) 
in-7G 
2 


where the factor U, called the “Overall Coefficient of Heat Transfer”, is a function of air ve- 
locity and, to a small degree, temperature and rate of condensation. 


It is common practice to show the results of air-cooled condenser tests by plotting this co- 
efficient U against the air velocity V, neglecting the effect of temperature and rate of con- 
densation. A much more useful way of summarizing test data is to give equations or charts 
for the values of the individual film coefficients which combine to form the overall coefficient. 
The relation between the coefficients is 


1 1 1 
7+ (2) 
U hg Rhy 
As shown in detail in Appendix 1, the heat dissipated by a wetted-surface, air-cooled con- 
denser can be computed in an analogous way by using the equation 


q=-u'a Sh - 4h (3) 


where the factor U', which may be called the “Overall Coefficient”, is again a function of the 
air velocity and, this time to an appreciable degree, the rate and temperature of condensation. 
An additional quantity which might affect the value of U'is the rate of circulationof water 
over the external surface, since this will determine the thickness, distribution and turbulence 
of the layer of water on the external surface and, in addition, may have a noticeable effect on 
the turbulence of the air stream. 


Thus, the results of tests on wetted-surface condensers can be shown by plotting curves for 
the coefficient U' against the air velocity V for a number of constant values of the other in- 
dependent variables or, exactly as in the case of dry-surface condensers, by giving equations 
or curves for the individual coefficients which this time can be combined to form U' by the 
equation 





Rhy, (4) 


where s and r are constants; q is a known function of the condensing temperature and rate; and 
hg and hy, are individual film coefficients of the same sort as occur in the equation for the 
overall coefficient for dry-surface condensers. For the sort of work with which this paper is 
concerned, these individual coefficients may be accurately enough represented by equations of 
the type 


hy =cv°* (5) 
hy, =D (6) 


In these equations the factors C and D will be constants, as in the case of dry-surface ap- 
paratus, unless the effect of changing the water rate is found to be important. In this case C 
and D should be expressed as functions of the rate of water circulation, thus including that 
variable in the equation for U'. 


Because there are so many independent variables it is not very practical to plot curves for 
U' directly from the test data or even to show the test results in this form. However, if the 
tests are separated into groups run at constant water rates (so as to avoid immediate con- 
sideration of this variable), it is fairly easy to determine the values of the individual coef- 
ficients by any of a number of methods. 
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In the present case, the test points (Table II) were not preselected for the application of the 
theory so that it was necessary to develop a method of analysis that did not require selected 
test points. The method is described in detail in Appendix 2 of this paper. Briefly, it consisted 
in assuming values for the individual film coefficient constants C and D and then using these 
assumed values to calculate values of the overall coefficient U' for each of the test operating 
conditions by using Eq (4). In general these values were not equal to the values computed 
directly from the test data by using Eq (3) in the form 


+ <5 ot (7) 
lh — 
Sk 
The differences between the values of U' computed by these two methods were said to be due 
to errors in the assumed values of the constants C and D. The amounts of these errors were 
then computed by the methods given in Appendix 2, and their average values were added to the 
originally assumed values of C and D to get new values for these constants, which procedure 
caused the equations to fit the test data to a satisfactory degree of precision. The originally 
assumed values of the constants were those which previous tests had yielded for dry-surface 
condensers. By referring to the discussion of the theory in Appendix 1, it will be seen that 
the values of the individual coefficients hg and hy should be about the same for wet- as for 
dry-surface condensers. 


In the present case, the application of these methods of analysis reduced the test data to the 
equations 


hy = 0.212 y’* (8) 
hy, = 200 (9) 


By inserting these specific values in the general equations, it is easy to synthesize a usably 
accurate description of the thermal performance of any proposed evaporative condenser. 


How well such a description fits the experimental condenser is shown in Figure 2 where 


values of heat dissipated computed by the above method are compared with the observed test 
values. 
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Figure 2 

When a number of designs are under consideration or when complete performance curves 
are to be made, a great deal of time and effort can be saved by plotting the equations with 
their test-determined constants in the form of nomographs or alignment charts. One such set 
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of charts based on the data of this paper is shown on Figures 3, 4 and 5. As an example of 
their use, suppose we have given q = 305,000 btu per hr, A = 410 sq ft, RA = 147 sq ft, t, = 
111.6 F, and V = 5,000 cfm, and it is desired to find what inlet air wet bulb temperature will 


be required. 
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Figure 3 


On Figure 3 draw a straight line from the given value, q/RA = 2070, on the left-hand scale 
to the given value, ty = 111.6, on the right-hand scale. At the intersection of this straight line 
with the middle scale read hy,/(s + ar) = 100. 


Multiply this by R, the given ratio of internal to external surface, and then on Figure 4 draw 
a straight line from Rhy,/(s + ar) = 35.8 on the right hand scale to the given quantity V = 400 
on the left hand scale. At the intersection of this straight line with the middle scale read 
U' = 16.55. 


Next, on Figure 5 draw a straight line from U' = 16.55 on the left-hand scale to the given 
quotient W/A = 12.2 on the right-hand scale. At the intersection of this line with the middle 
scale read q/(4i,W) = 1.17. 


Ai is equal to the quotient of the given q/W divided by this quantity, 1.17, and the value of 
i, the inlet air enthalpy, can be obtained by subtracting this quotient, Ai, = 52.1, from 88.1 
the enthalpy of saturated air at the refrigerant temperature. Having thus determined that 
i, = 36.0, the corresponding wet bulb temperature of the inlet air, tj; = 80 F, can be read from 
any psychrometric chart or table. 


In case the only interest is in exploring design possibilities, the labor can be greatly de- 
creased by using constant values for the refrigerant condensing temperature and inlet air wet 
bulb temperature. This, in effect, sets a general design standard for the temperature differ- 
ence. When this is done the remaining variables can all be shown ona single nomograph so 
that all the relations between the major variables for a complete design can be determined 
simply by drawing one straight line. 
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Figure 5 


Such a nomograph, based on the data of the present tests, is shown in Figure 6. This nom- 
ograph shows the relation between heat dissipated, External Surface, Internal Surface, Air 
Volume, and Air Velocity for fixed values of Inlet Air Wet-Bulb Temperature and Refrigerant 


Temperature. 
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Figure 6 


A typical design study is shown on Figures 7 to 11. These exhibit the possibilities in the 
design of an evaporative condenser for a twenty-ton unit. The design performance was set at 
153,000 btu/hr at a condensing temperature of 111.6 F and an inlet air wet-bulb temperature 
of 80 F. 
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Figure 7. Design Study of a Freon-12 Evaporative Conden- 
ser Showing the Balance Between Surface and Air Volume 
for the Design Conditions Condensing Pressure - - 154 psia; 
Inlet Air Wet-bulb - - + 80F; Heat Dissipated - 153,000 
btu/hr. 
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Figure 7 was plotted from the nomograph of Figure 5. From a study of this figure it is easy 
to make a preliminary choice of the quantity of air to be used through the condenser, because 
an increase in the product U' A is equivalent to an increase in cost. With this in mind, it will 
be apparent that the proper air volume will be located near the knee of the curve. A reason- 
able preliminary choice would be 2,500 CFM. 


Figure 8 shows the relation between internal and external surface areas for a number of air 
velocities, all at the previously selected air volume of 2,500CFM. These curves were plotted 
from the nomograph of Figure 6. 





hateroal Area 


Figure 8. Design Study of a Freon-l2 Evaporative Condenser 
Showing the Balance Between Internal and External Surface for 


the Design Conditions: Condensing Pressure - + 154 psia; 
Inlet Air Wet-bulb~- - + 80F; Heat Dissipated - 153,000 btu/hr; 
Air Volume + + + + + 2,500 cfm, 


By making a tentative selection of materials for the tubes and fins the curves of Figure 8 
can easily be converted into curves of weight against external surface. The results of sucha 
conversion are shown in Figure 9, which is based on the use of copper tubes of 0.035-inch wall 
thickness and copper fins 0.022 inch thick. From this figure, it is at once evident that for any 


given air velocity there will be one particular design that will result in the most effective use 
of material. 


Figure 10 was plotted directly from Figure 9 to show how the minimum weights vary with 
air velocity. Each point on the curve of this figure represents the lightest condenser, built of 
the specified materials, that can be designed to give the required performance at the indicated 
air supply. While no minimum appears on these curves, it should be bourne in mind that the 
cost of fan and motor for supplying the air increases as the air velocity is increased so that 
if the weights were converted to costs and added to the cost of fan and motor the resulting 
curve of cost against air velocity would certainly show a minimum. 


Figure 11 shows how the various physical dimensions of the condenser of minimum weight 
vary as the air velocity is changed, all at an air supply of 2500 CFM. 
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Figure 9. Design Study of a Freon-l2 Evaporative 
Condenser Showing the Balance Between Surface, Air 
Velocity and Weight of Active Material for the Design 


Conditions: Condensing pressure - + 154 psia; Inlet 
Air Wet-bulb - - - 80F; Heat Dissipated - 153,000 
btu/hr; Air Volume - - - - - 2,500 cfm; Tube Mate- 


rial 0.035 inch Copper; Fin Material 0.022 inch Copper. 
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Figure 10 


APPENDIX 1 
Derivation of the Equations for Heat Flow 


Consider the simple wetted-surface heat transfer apparatus shown in Figure 12. It consists 
of a rectangular, insulated duct of constant width y and length x divided lengthwise by a thin 
metal wall. The space below this wall is kept full of a pure saturated vapor which is con- 
densing on the wall, while the portion of the duct above the dividing wall serves as a channel 
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Outside Wall of Duct (made of insulating material) 







Stream of Moist Air Above Wall 


Thin Layer of Liquid Water 












Film of Condensate 


mdensing Refrigerant Vapor Below Wall 


if Duct Section A-A 


Simplified Wetted-surface Condenser Used in Developing the Theory of Heat Flow. 


for a stream of moist air. The air-side of the dividing wallis kept covered with a thin film of 
liquid water. The refrigerant or vapor side evidently will be covered by a film of condensate 
from the saturated vapor. 


When the moist air is at or near normal atmospheric pressure and the temperature of the 
dividing wall is less than 212 F, it is obvious that the hydraulic pressure exerted on the water 
film by the air will prevent the formation of vapor bubbles in the water, so that all evaporation 
must take place at the surface of the water film. Therefore, all the heat released at the surface 
of the condensate film by the condensation of the refrigerant must pass through the condensate 
film, dividing wall and water film, in series, to the surface of the liquid water. Here, part of 
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the heat will be used to vaporize the water while the remainder passes on through the “stagnant 
air film” to the main body of the air stream. The relative rates of flow of heat and diffusion 
of water vapor through this air film will determine what part of the heat is used to evaporate 
the water, thus increasing the moisture content or humidity of the air stream, and what part 
passes through to the air stream directly, thus increasing its temperature. 


The rate at which any small element ydx of the surface of the water is gaining heat from the 
refrigerant is 

dq = (t, - tw)hyydx (1-1) 
in which the condensate film, metal wall and water film are considered a composite wall, that 
is, 


— See 4¢ = + == (1-2) 


The rate at which this element of surface is losing heat to the air stream directly is 
dq, = (ty - t)hpydx (1-3) 


while the rate at which water vapor is diffusing from the water surface through the air film to 
the main air stream is 


dw, = (Hy, - H)Kydx (1-4) 


The heat absorbed by the evaporation of this quantity of water is 
dq, = rdWy (1-5) 
so that the rate at which the element of surface of the water is losing heat by evaporation is 
dq, = r(H,, - H)Kydx (1-6) 
In the steady state, when conditions at all points are independent of the time, it is obvious 
that heat cannot be stored anywhere in the apparatus. Hence, the rate at which the element of 


water surface is losing heat to the air stream by conduction and evaporation must be equal to 
the rate at which it is gaining heat from the condensing refrigerant so that 


dq = dq, + dq, = [(ty - t)hy + r(Hy - H)K] ydx (1-7) 


As a matter of convenience write 


K = hy/B (1-8) 
Then after substituting this value in Eq (1-7) and rearranging terms we have 
dq = F [( Btw + rHy) - (Bt + rH)) ydx (1-9) 
But from Eq (i-1) 
4a (1-10) 


and if we write 


So that after using Eqs (1-10) and (1-12) to eliminate the conditions at the water surface from 
Eq (1-9), we obtain after grouping terms 


dq = U' Pydx (1-13) 


where as a matter of convenience we have written 


U' hy hy 


1_ 8, B+ar (1-14) 


and 
P = (Bt, + rH,) - (Bt + rH) (1-15) 


Integrating Eq(1-13) between the limits yx = 0 and yx = A and using the Mean Value Theorem, 
we obtain 
qe] U' P ydx = Umean { P ydx = Umean A Pmean (1-16) 


In the case of air and water vapor mixtures ) is very nearly equal to s so that, for work of 
this sort, 4i can be substituted for P. Also, experience has shown that the logarithmic mean 
is a convenient and practically useful approximation to Pmean- When these approximations 


are used, it is easy to compute corresponding values of U}mean directly from the test data by 
the equation 


Umean * big ‘ Ai, (1-17) 





— = (1-18) 
U's Rhy 


where a is to be evaluated at the mean value 


= te - Rhr (1-19) 


The constant R has been introduced as a multiplier of hy; in Eqs (1-18) and (1-19) to permit 
expressing this coefficient per unit area of refrigerant side surface when dealing with finned 
tubing. In this case the areas of the air side and refrigerant side of the wall would not be 
equal, as they were assumed to be in the preceding discussion. 


An estimate of the error involved in the use of the logarithmic mean can be obtained in the 
following manner. 


By equating the heat absorbed by the air in passing over an element of the surface to the heat 
transfered to it from the surface we obtain 


di 
Gap dP = U P dx (1-20) 
so that, taking di/dP = 1, 
—_ 
p-pPexp{ Dax (1-21) 
0 G 
Let 
(*u; 
P' = P exp) —'dx (1-22) 
0 G 
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The logarithmic mean is 


a j 
log mean P otra 3 } P, exp Gs P dx (1-23) 
inp i 
But by integration of Eq (1-21), 
a 
Pp. U' } A 
ets) & . YVmean e 
n P, ) G dx aarenaee (1-24) 
Therefore 
a a 
1 lx U' 
log mean P =— } P.exp |= = ax) dx (1-25) 
A . \A } G 
P oS ft om (2 dx) dx (1-26) 
mean * x } ,exp } G } J 
and a - \ 
' 1 U' 
Pmean = A } P, exp({ G ax} dx (1-27) 


From which since U' continuously increases as x increases, it is obvious that 
Pmean <“Pmean <!0g mean P (1-28) 


so that the error incurred by using the logarithmic mean will be less than that given by the 
equation 


Pmean = (1 - E) log mean P (1-29) 


where 


(By -1 


s .< 1. (1-30) 


e) 


1 


with 


(1-30) 


APPENDIX 2 
Evaluation of the Individual Coefficients 


Assume values for C and D in Eqs (5) and (6). Then Eq (4) can be written 

1 s f1 [1] [1 | 

Ht. ec sane IC }o r (Rhy J - r(a), LRhy J, ole 
where the subscripts , are used to denote assumed values of C and hy, and the corresponding 
calculated value of U' as distinguished from the test value of U' computed from the test data 
by means of Eq (7). The difference between the calculated value and the test value of U' for 


each test can be regarded as due to errors AC and jhy, in the assumed values of these con- 
Stants. That is 





1 ee: 1 1 | 
U'~ LY o= ys Ag+ (s+ar), ORR + Tl Rh | Ad (2-2) 








For any constant value of t, and over a small range of values of t, - ty we may write as a 
close approximation 


40 = aA(ty - ty) (2-3) 


where a is a constant. This is equivalent to using a quadratic equation, H = a, + a,t + a,t, to 
fit a short section of the t-H curve. The constant a, in the quadratic equation would then be 
equal to the constant a in Eq (2-3). 


Since 
| 
tr - ty = Rhy. (2-4) 
and since q/A will have a fixed value for each test, we can write, 
1 
Sa=a(q/A)A Rhy, (2-5) 
With this substitution Eq (2-2) becomes 
a yy q/Aj| , 1 
A pz yaAlsct (s+ ar), + ar [a4] 4 Rh, (2-6) 
so that each test point supplies an equation for the variables / (1,4) and A(1/Rhz) of the form 
1 1 
A.A G+ Ad Rhy, ~ A; (2-7) 


Where the constants A,, A, and A, will have different values for the various test points. Their 
values can be computed from the test data by means of the equations, 





A, =a33 (2-8) 

A, = (s+ ar), + ar tical (2-9) 
e- 2 ph) 

saga tel, — 


Thus, any two tests will furnish a pair of simultaneous equations from which the values of 
4(1/Rhz,) and A(1/C) can be calculated. The results of more than two tests can be averaged 
by the method of least squares. In this case the pair of simultaneous equations, called the 
“Normal Equations”, are 


ae Za? + re ZA,A, = EAA, (2-11) 
: 1 sa? 2-12 
4¢ *A A+ Spy DA, = ZA,A, (2-12) 


in which each term of the summation includes data from only one test. 


If it is desired to assign different weights, W,, W, etc., to the various tests, this can be 
done by adding the data from each test W times in making the summations, that is, the sum- 
mation terms will be EWA‘, EWA,A, etc., instead of EAi, ZA,A, etc. For example, in the 
case of the tests at 2500 - 2600 lbs per hr water rate, 10 of the 13 tests were run at an air 
velocity of about 300 fpm. In order to avoid overemphasizing the 300 fpm velocity, each of 
the 10 tests at this velocity was assigned a weight of unity, while each of the remaining tests, 
at other air velocities, was given a weight of 10. 


The assumed values of C and D used in applying this method of analysis to the data of Table 2 
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were based on the results of previous tests on dehumidifiers and dry-surface condensere. 
The values were 


(C), = 0.178 and (D), = 200 (2-13) 

that is 
(hg), = 0.178 V°"* and (hy,}, = 200 (2-14) 

The results of the analysis of the data of Table 2 are as follows: 
Test Group I - Water rate 3700 - 3800 lb per hr 
All tests given equal weight, 

hg = 0.237 V°* and hy = 212 (2-15) 
Test Group II - Water rate 2500 - 2600 lb per hr 
A - All tests given equal weight, 


hy = 0.251 V°" and hy = 178 (2-16) 


u 
B - Each of the 10 tests at 300 fpm air velocity given a weight equal to 0.1 that of each of 
the three tests at other air velocities. As previously mentioned, the purpose of this was to 
avoid overemphasizing the 300 fpm velocity. 


hg = 0.212 V° and hy, = 190 (2-17) 


After considering these results of the analysis it was concluded that over the range 2500- 
3800 lb per hr the effect of a change in the water rate could be neglectedand the test data could 
be represented by the equations 


hy = 0.212 v°* and hy, = 200 (2-18) 


g= 
Figure 2 compares the values of U' calculated from these coefficients with the values ob- 
tained directly from the test data. 


APPENDIX 3 
List of Symbols 


A = Total area of heat transfer surface exposed to air - sq ft 

A,, A, and A, = Constants defined by Eqs (2-8), (2-9) and (2-10) 

a = Proportionality factor defined by Eq (2-3) 

C = Constant defined by Eq (5) 

D = Constant defined by Eq (6) 

G = Mass rate of flow of dry air over heat transfer surface-lb of dry air per hr 

H = Average humidity of moist air over cross-section of stream at x - lb of water vapor 
per lb of dry air. 


Hr, = Humidity of saturated air at condensing temperature “tr” 

Hy = Humidity of saturated air at temperature “ty” 

H, and H, = Humidity of air entering and leaving heat transfer apparatus 

he = Coefficient of heat transfer from condensing refrigerant vapor to heat transfer surface 
- btu/(hr) (sq ft) (F) 

hg = Coefficient of heat transfer from heat transfer surface to air - btu/(hr) (sq ft) (F) 

ht, = Coefficient of heat transfer from condensing vapor to air side surface of water film 


- btu/(hr) (sq ft of vapor side surface) (F) 
i = Average enthalpyof moist air over cross-section of stream at x - btu per lb of dry air 
ir = Enthalpy of saturated air at condensing temperature “ty 
i, and i, = Average enthalpy of air entering and leaving heat transfer apparatus 
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tr = 
tw = 


Coefficient of diffusion of water vapor from water surface to air - (lb of water)/(hr) 
(sq ft) (unit humidity difference between saturated air at temperature of water surface 
and air) 


Thermal conductivity of dividing wall - btu/(hr) (sq ft) (F/ft) 

Thermal conductivity of liquid water 

Rate of flow of heat from condensing refrigerant vapor to air stream - btu per hr. 
Ratio of area of heat transfer surface exposed to refrigerant to area of heat transfer 
surface exposed to air 

Latent heat of evaporation of water - btu per lb 

Humid heat of moist air (i.e.,the specific heat per lb of dry air) - btu/(F)(lb of dry air) 
Average temperature of moist air over cross-section of stream at x - F 
Temperature of condensing saturated refrigerant vapor - F 

Temperature of water surface in contact with air - F 


t, and t, = Average temperature of air entering and leaving heat transfer apparatus 


U = Overall coefficient of heat transfer for dry-surface air-cooled condensers - btu/ (hr) 
(sq ft of air side surface) (F) 
U' = Overall coefficient for wetted-surface air-cooled condensers -1lb/(hr)(sq ft of air side 
surface) 
V = Average linear velocity of air stream - ft per minute 
Wy = Mass rate of diffusion of water vapor from water surface to air - lb per hr 
a = Proportionality factor defined by Eq (1-11) 
sp = Thickness of dividing wall - ft 
dw = Thickness of water film - ft 
44, = ip - 1h 
Ai, = ir - iz 
4t, = tr-t, 
bt = t.-t, 
TABLE 1 
Details of Test Condenser 
Core Size Tube Details 

Height 2.14 ft. Material Copper (tinned) 
Width 1.00 ft. Outside diameter 0.500 in. 
Depth 1.21 ft. Wall thickness 0.028 in. 
Face Area 2.14 sq. ft. No. of rows 9 

No. per row 16 

Spacing 1.56 in.x 1.56 in. 

Surface Area Fin Details 

Tubes (External ) 18.8 sq. ft. Material Copper (tinned) 
Fins 65.3 sq. ft. Height 2.14 ft. 
Total External 84.1 sq. ft. Width 1.21 ft. 
Internal 16.7 sq. ft. Thickness 0.022 in. 
Ratio 1 to 5.04 Spacing 0.75 in. 





Does not include return bends. 
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Table II - Test Data 






















































































; CO Heat Released 
a Air Water B.t.u./hr. 
| ; ; Temp. Note 1 | 
| Eh. 1. a a Rate Cond. | Calculated | Calculated 
Inlet | Outlet 
| Velocity Spray of Pr. from from 
Test fpm Dry | Wet| Dry | Wet |Nozzle |Circulation |Lb/sq.in.| Refr. Air % 
| Nc | i Bulb |Bulb| Bulb; Bulb; OF Lb/Hr Absolute | Data Data Differ. 
3 | 303 100.0} 80.2} 99.1) 94.9] 97.0 2570 162.9 | 50000 52800 5.6 
4 | 305 | 90.0/80.1) 96.3] 94.5| 96.2 2570 160.4 | 49600 51450 3.7 
6 305 110.0/ 94.0 | 108.0/103.0 | 106.3 2570 | 180.4 46500 42500 -8.5 
7 | 303 | 90.0/73.5| 93.5] 90.5] 92.3 2570 =| 154.4 | 49300 53400 8.2 
8 | 305 | 90.0/83.3| 97.2] 96.6) 98.2 2570 164.4 | 48400 51050 3.7 
9 305 |100.0192.2 }103.5|102.9 | 104.1 | 2570 174.9 46950 49550 5.7 
10 | 305 |100.0|86.0|100.5) 98.8 / 100.2 | 2570 | 167.9 48100 | 52100 8.3 
11 307 | 110.0/80.2 |100.8 94.9| 96.6 2570 | 162.0 49100 | 53600 9.1 
12 306 |100.0)75.4 95.7| 91.7} 94.1 | 2570 | 155.9 49000 | 54150 10.6 
1a | 317 | 97.0/76.7| 96.6| 92.5| 94.6 | 2570 | 156.9 49450 | 55550 | 12.3 
2A | 309 97.0|76.4 |106.8) 94.2] 98.3 | 780 | 196.9 45900 | 61500 | 34.0 
4A | 305 97.0/76.8 | 94.5) 91.5| 95.2 3800 153.5 50150 49450 | -1.5 
ic | 319 100.0/80.9 | 97.0) 95.0} 98.0 3780 158.0 49450 | 54450 | 10.0 
2C 319 100.0/80.9 97.3) 95.0 | 98.4 2920 162.5 | 48700 54450 | 11.7 
3C | 319 |100.0/80.9 |102.6| 97.0| 102.4 | 1680 182.5 44450 63800 | 43.6 
6C 315 100.0|83.2 100.6} 96.8 100.8 2125 171.5 47500 93500 | 12.8 
7 | 319 100.0 /80.5 |108.1) 94.5 | 105.2 692 } 215.5 41800 53450 | 28.0 
1D 319 100.0)80.0 96.0| 93.2} 96.5 | 3770 | 157.0 48750 49150 | 1.0 
2D 131 100.0|80.0 |105.9/105.0/107.1 | 3780 | 178.5 45950 44700 -2.7 
3D | 194 100.0 /80.0 |101.8|100.0 | 102.6 | 3770 169.6 46000 | 48600 5.9 
4D | 466 100.0 |80.4 | 95.4] 90.0} 93.9 | 3660 | 151.0 50400 | 49300 -2.1 
1E | 128 |100.0/79.7 |106.8|105.9|108.8 | 2570 187.0 44800 | 46500 | 3.8 
| 2E | 322 100.0|79.7 98.8) 93.8 | 98.0 | 2535 163.0 48700 | 52900 8.4 
| 3E | 481 100.0 /80.0 | 96.8 | 90.8 | 94.9 | 2525 156.6 ‘| 49300 | 58500 18.8 
Note 1 - The values of Q used in the analysis were taken from this column. 
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